Cluster packing geometry for Al-based F-type icosahedral alloys 
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This paper presents a new highly stable periodic approximant to the Al-based F-type icosahedral 
quasicrystals, i-Al-Pd-TM (TM=transition metals). The structure of this intermetallic Al-Pd-Cr- 
Fe compound is determined ab initio using single-crystal X-ray diffraction, where the space group 
is identified to be Pa3 and the lattice constant 40. 5A. The structure is well described as a dense 
packing of clusters of two kinds, which are known in the literature as the pseudo-Mackay type and 
the Bergman type clusters. The clusters are centered at the vertices of a canonical cell tiling, in 
which the parity of each vertex determines the kind of the associated cluster. Adjacent clusters can 
be markedly interpenetrated, while the structure requires no glue atoms to fill in the gaps between 
the clusters. It is shown that the crystal can be designated as a 2 x 2 x 2 superstructure of the 
ordinary cubic 3/2 rational approximant. The superlattice ordering is shown to be of a different 
kind from the P-type superlattice ordering previously reported in i-Al-Pd-Mn. The present results 
will greatly improve the understanding of atomic structures of F-type icosahedral quasicrystals and 
their approximants. 

PACS numbers: 61.44.Br, 61.66.Dk, 61. 50. Ah, 61.05.cp 



I. INTRODUCTION 

The Al based F-type icosahedral quasicrystals, such 
as i-Alg5Cu2oFei5 and «-Al7oPd2oTMio (TM=transition 
metal, e.g., Mn, Re), pose a number of questions as to the 
detailed atomic arrangements. Researchers have tried to 
extract as much structural information as possible from 
materials of this kind through state-of-the-art techniques 
for structure analysis using X-ray, electron or neutron 
diffraction. These studies indicated that there existed a 
few different kinds of clusters as the basic building units 
that construct the quasicrystals. 

A direct measurement of a quasicrystalline sample, 
however, entails a well known disadvantage in terms of 
structure determination. Conventional methods devel- 
oped for analyzing crystal structures are inherently inap- 
plicable to quasicrystals due to the absence of periodicity. 
However, these methods could still provide invaluable in- 
formation regarding the constituent clusters as well as 
their local packing in quasicrystals through analyzing 
their rational approximants. In particular, the latter ap- 
proach has so far been taken successfully in elucidating 
local characteristics of P-type icosahedral quasicrystals, 
such as i-Al 73 Mn 2 iSiei' 2 - and i-Cd 5 .7YfA 

Unfortunately, for the case of F-type icosahedral qua- 
sicrystals, a stable approximant phase has scarcely been 
reported. Hence, local characteristics of the structure re- 
mained to be uncertain to a large extent. It appears 
as if Al-based F-type quasicrystals are so stable that 
no subordinate approximant phase could be obtained 
via a slight change of the composition. To date, only 
two possibly related cubic approximants have been pre- 
viously reported in a rapidly quenched alloy with a nom- 
inal composition of Alg9Pd2oMngSi3 after annealing at 
above 1000°C 4 . By way of single crystal X-ray diffrac- 
tion, these approximants were determined to be cubic 



1/1 and 2/1 approximants to the F-type icosahedral qua- 
sicrystal based on their lattice constants of 12.28A and 
20.2lA— The structural information of the two refined 
structures has been utilized as the basis for building a 
six-dimensional structure model of the relevant F-type 
icosahedral quasicrystal 7 . However, the reliability of the 
two structures has remained controversial because their 
constituent clusters were totally different despite their 
compositional similarity. 

In i-Al7oPd2oTMio and its close associates, Al, Pd and 
TM are likely to maintain their distinct roles in the struc- 
ture. While no favorable approximant to the icosahedral 
quasicrystal has been reported for TM = Mn or Re, an 
ample room remains for the choice of TM elements. In 
particular, a fine adjustment of the electron concentra- 
tion (e.g., e/a ratio) could be achieved by mixing two 
transition metal elements for TM. In this report, a syn- 
thesis of a new stable approximant is accomplished by 
blending Cr and Fe, the two immediate neighbors to Mn 
in the periodic table, for TM. An optimal ratio between 
Cr and Fe is searched experimentally The stability of 
the new phase is such that fine single crystals may grow 
through a simple slow cooling method. Single crystal X- 
ray diffraction is then performed to analyze the structure. 
The first part of this paper is devoted to a presentation 
of the synthesis and the crystal structure analysis. 

The present structure analysis offers a source of infor- 
mation that is essential in understanding the local atomic 
arrangements in the F-type icosahedral quasicrystals. It 
turns out that the structure is composed of two kinds 
of clusters, called the pseudo-Mackay type and Bergman 
type clusters. The clusters are interconnected in such a 
way that significant interpenetrations are allowed, while 
no glue atom between them is required. Moreover, the 
centers of the clusters are given as the nodes (or vertices) 
of a three-dimensional tiling with four kinds of polyhedra 
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called the canonical cellsS; this finding is of fundamental 
importance as it allows a systematic description of the 
basic skeleton of this as well as other related compounds. 
The main body of the work comprises a full account of 
the crystal structure and related discussions. 

The structural description can be generalized in a 
straightforward manner, leading to the proposal of a 
number of hypothetical structures which could form as 
real approximants. Although a full description of the F- 
type icosahedral quasicrystals awaits future endeavors, 
our understanding on the local atomic arrangements in 
the F-type icosahedral quasicrystals and their approx- 
imants can now be revised significantly based on firm 
experimental evidences. This yet opens a way to further 
attempt to synthesize variants of stable approximants as 
well as to gain insights on the superlattice ordering phe- 
nomenon which was reported previously in the F-type 
icosahedral quasicrystal i-Al7oPd2oTMi<j£. 

This paper is organized as follows. Section fJTT] illus- 
trates how a simple compositional search culminates in 
the discovery of a new approximant phase. The initial 
samples were polycrystalline, yet powder X-ray diffrac- 
tion as well as electron diffraction proves the presence of 
a cubic approximant with a large unit cell. An electron 
micro-probe analyzer (EPMA) is used to evaluate the 
fine chemical composition of the approximant. Then an 
optimized synthesis is performed, resulting in the growth 
of single crystals exceeding 100 fim in diameter. Single 
crystal X-ray diffraction and the structure analysis are 
performed in Section fJTTTJ The crystal structure is de- 
scribed in detail in Section ^IVI in a constructive manner: 
After introducing the geometrical templates of the two 
kinds of clusters, the packing geometry of the clusters is 
described in detail. It is shown that the global arrange- 
ment of the clusters is determined based on a canoni- 
cal cell tiling with F-type ordering. In an early part of 
Section SjV] discussions are given of the relationship be- 
tween a superlattice ordering of the present approximant 
and that in the quasicrystal. A likelihood of anti-phase 
boundaries as a possible source of disorder in the present 
material is also discussed. Then the rest of this section 
examines the possibility of introducing a tiling model for 
the atomic arrangement. Section i jVIl is devoted to con- 
cluding remarks. 



II. SYNTHESIS AND CHARACTERIZATION 

Our synthesis of the new approximant was performed 
in two steps. Firstly, a search for a possible approximant 
phase was made by varying the ratio between Cr and Fe. 
A few characterization techniques were used to locate the 
composition of the approximant phase as well as to ex- 
tract the basic crystallographic information. Secondly, 
starting from the right composition a new synthesis was 
performed with a slow cooling. High quality single crys- 
tals of the approximant were successfully obtained in this 
way. These steps are described in the following. 



FIG. 1. Powder X-ray diffraction diagrams for 

Al7oPd2oCr I Feio-a; (x = 3, 5 or 7). For x — 3 and 5, 
the peaks marked with filled circles are associated with an 
approximant, while for x = 7 those marked with filled squares 
are associated with a quasicrystal. Note that unassigned 
peaks associated with impurity phases become prominent as 
x is increased 
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The starting materials were Al (Kojundo Chemical 
Lab.; purity, 99.9%), Pd (Tanaka Kikinzoku; purity, 
99.95%), Cr (Furuuchi Chemical; purity, 99.9%) and Fc 
(Nilaco; purity, 99.5%). An alloy ingot with a nomi- 
nal composition of Al7 Pd2oCr a ;Feio- I [x = 3, 5 or 7) 
was prepared by the arc-melting method under argon at- 
mosphere. After sufficient homogenization was achieved, 
the as-solidified ingot was fragmented, put into a Tam- 
mann crucible made of AI2O3 and sealed altogether into 
a quartz tube with pure argon gas of about 0.08 MPa. 
In an electric furnace, the sample was annealed at 850°C 
for 48 hours and cooled swiftly (~ 1 hour) down to room 
temperature. 

In Figure [T] powder X-ray diffraction patterns taken 
from the samples x = 3, 5 and 7 are shown, where 
the characteristic X-ray of Cu Ka (A = 1.543A) is used 
with Bragg-Brcntano diffraction geometry (Mac Science, 
diffractometer M03XHF 22 ). Slightly above the strongest 
Bragg peak at 29 ~ 44°^ a clear sub-peak is observed for 
the two samples x = 3 and 5. The latter peak feature 
is usually associated with the cubic 2/1 approximant to 
an icosahedral quasicrystal with the Miller indices being 
hkl=10 0. The corresponding lattice constant is calcu- 
lated to be c.a. 20. 3A. The sample x — 3 appears to con- 
tain the largest amount of the approximant phase, while 
the suppression of the sub-peak as well as the growth of 
an extra feature at 29 ~ 43° clearly indicates the preva- 
lence of impurity phase(s) as x is increased. Hence, it is 
expected that Al7oPd 2 oCr 3 Fe7 is the closest composition 
to that of the approximant phase. 

Back-scattered electron images taken from the three 
samples using a scanning electron microscope (SEM) 
are shown in Figure [2] The SEM apparatus used is 
JXA-8621MX (JEOL), in which an electron-probe micro- 
analyzer (EPMA) is implemented for chemical compo- 
sition analysis. A few domains of different gray lev- 
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FIG. 2. SEM back-scattered electron images for (a) 
x=3, (b) 5 and (c) 7. The estimated composi- 
tions are: (a) Al69.iPd22.oCr2.iFe6.s (light gray) and 
Al 6 8.4Pd2o.9Cri.iFe 9 .5 (dark gray), (b) Al 6 9.7Pd22.5Cr 2 .3Fe5.4 
(light gray) and Al7 2 .sPdii. 8 Cri 1 . 6 Fe4.o (dark gray) and (c) 
Al69.5Pd23.5Cr3.6Fe3. 4 (light gray) and Al 3 Pd 2 (white). 




els are observed in each sample (voids are shown as 
black regions), indicating the existence of a few alloy 
phases with different compositions. It is observed that 
a light gray region in Figure [2ja) represents the majority 
phase for x = 3 with the estimated composition being 
Alg9.1Pd22.0Cr2.1Fe6. 8- It is reasonable to assume that 
this region corresponds to the approximant, which con- 
stitute the main body of the sample. For x — 5, a light 
gray region in Figure [3Jb) has an estimated composition 
of Al69.7Pd22.5Cr2.3Fe5. 4, and one can still associate it 
with the approximant for the compositional difference is 
insubstantial. However, all the phases observed for x = 7 
are deviated significantly from the approximant in terms 
of chemical composition, implying no relevant approxi- 
mant phase exists in the sample. 

Firmer evidence of the approximant phase was pro- 



FIG. 3. TEM diffraction patterns taken from the 
Al7oPd2oCr3Fe7 sample along (a) the two-fold axis (|| < 
100 >), (b) the three-fold axis (|| < 111 >) and (c) one of 
the pseudo five-fold axes. 
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vided by selected-area electron diffraction. Here two 
transmission electron microscopes (TEM), JEM-2000EX 
(JEOL) and JEM-2000EXII (JEOL), operating at 200 
kV were used. Electron diffraction patterns taken from 
the Al 70 Pd2oCr 3 Fe7 sample (Figure [3]) clearly indicate 
the existence of 2-fold, 3-fold and pseudo 5-fold axes, 
where strong Bragg reflections are arranged in a similar 
way to the case of the quasicrystaliS. In Figure OUa) and 
(b), periodic arrays of spots can be indexed as a cubic 
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FIG. 4. SEM micrograph taken from the single crys- 
talline sample obtained with a nominal composition of 
Al69.4Pd22.3Cr2.2Fe6.!. 
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crystal with a lattice constant of about 40 A. Note that, 
however, this is twice as large as the one estimated from 
the above powder X-ray diffraction pattern by assum- 
ing a conventional 2/1 approximant. Moreover, extinc- 
tions are observed at Okl with fc=odd and at 00/ with 
Z=odd, indicating that the crystal has a non-symmorphic 
space group. Therefore, the pre-existing model of a 2/1 
approximant^iil clearly fails to describe the present ma- 
terial. This motivates us to carry out a thorough investi- 
gation of the crystal geometry based on an ab initio struc- 
ture determination using single crystal X-ray diffraction. 

In order to obtain a single crystal to be used for the 
structure analysis, a new sample was prepared in the fol- 
lowing way. A refined composition given by the arith- 
metic mean, Al69.4Pd22.3Cr2.2Fe6.ij of the two composi- 
tions evaluated from the samples x — 3 and 5 was used for 
the starting composition. The preparation steps were the 
same as those described above except that at the anneal- 
ing step within the furnace the sample was fully melted at 
1160°C for one hour and then cooled down slowly to room 
temperature with a cooling rate of 10°C/hour. Figure 
|4] shows an SEM micrograph taken from the new sam- 
ple using JCM-5100 (JEOL), where faceted crystalline 
grains roughly in the shape of a truncated octahedron 
are wrapped partially with debris (impurities). It con- 
firms that single crystals a few hundred /jm's in diame- 
ter were successfully grown, which also implies that the 
approximant formed congruently from the melt. 



III. AB INITIO STRUCTURE 
DETERMINATION 

A. Single crystal X-ray diffraction 

A single crystal was taken from the new sample with 
the estimated dimensions of 0.272mm x 0.178mm x 



0.320mm. The diffraction experiments were performed 
using a Bruker SMART APEX diffractometer, mounting 
CCD area detector, with Mo Ka radiation and graphite 
monochromator (wave length = 0.71073A). Indexing and 
empirical absorption correction were performed using the 
Bruker software package (SMART for WNT/2000 5.625 
Bruker AXS and SAINT 6.45 SADABS). In collecting 
the reflection intensities, the crystal was identified as a 
primitive cubic crystal with a lattice constant of 40.54A. 

Inspecting the original reflection data, it was found 
that a set of symmetrically equivalent reflections often 
included one or a few members showing exceptional de- 
viation in the intensity. This is due to the dynamical ef- 
fect (or multiple scattering) by which an extra intensity 
is added to the 'kinematical' Bragg intensity. Although 
the extra intensities are weak in contrast to those of the 
strongest Bragg reflections, they could still be harmful 
for weak reflections; note that the structural information 
carried by weak reflections plays a crucial role in deter- 
mining long-range characteristics of a complex structure 
such as a quasicrystal and an approximant with a large 
unit cell. 

In principle, a reflection intensity affected by multi- 
ple scattering could be effectively screened out from the 
data based on the intensity distribution among the set of 
equivalent hkl entries. In a common structure refinement 
software (e.g., JANA2006), the procedure can be per- 
formed when equivalent reflections are averaged. How- 
ever, if the number of equivalent hkl entries is small, the 
statistics is simply not enough to discern the ill-entry. 
Hence, we took a precaution to remove those hkl entries 
which have less than six symmetrically equivalent asso- 
ciates from the original dataset in order to avoid the pos- 
sible source of error. The preprocessed dataset was used 
for the structure analysis as described in the following 
subsection. 



B. Structure analysis 

Our structure analysis relied entirely on the JANA2006 
software packageA^. The Laue group and the space group 
was determined unambiguously to be m3 (Th) and Pa3, 
which turned out to be the only non-symmorphic space 
group with minor contradictions with the data. Out of 
the whole 95213 observed reflections satisfying / > 3<t(I), 
302 symmetrically extinct reflections were included. The 
symmetry averaging was performed while at the same 
time those data entries showing significant deviation from 
the averages of their equivalent reflections ( | / — I av | > 
10a(Ia V )) were eliminated from the data. The averaged 
data containing 8454 independent reflections satisfying 
I > 3tr(7) were obtained with R int = 7.91%. 

The initial structure model was generated with the SU- 
PERFLIP program, which is an implementation of the 
charge flipping algorithm in JANA2006. Then Fourier 
synthesis and least squares fitting were iterated. After 
several iterations, the automated peak search turned un- 
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TABLE I. Crystallographic data 



Formula 
Molar mass 

Temp, of data collection 

Space group 

Lattice constant, ai at 

Cell volume, O 

Z 

Calculated density 

Absorption coefficient 

Range of 26 

Independent reflections 

Obs. reflections (/ > 3<j(I)) 

-Rmt (obs/all) 

Num. of parameters 

R(F) 

R W {F) 

S 

III :.\ * 

Ap 

mm 



A/e.s.d. 



AI72.515 Pd22.498 Cr4.928 

5045.2 g/mol 
room temp. 
Pa3 (No. 205) 
40.5405 A 
66629.6 A 3 
40 

5.028 g/cm 3 
9.294 mm -1 

I. 74 ~ 53.2 
21526 

8454 

7.91/13.37 
783 

II. 96 
12.13 
4.59 

8.01, -5.38 e/A 3 



Fc 7 



0.0002 



successful in locating new atoms. Then we carefully in- 
spected the tentative structure, finding incomplete icosa- 
hedral clusters orderly packed in the unit cell. We 
identified two kinds of clusters, which in the literature 
are known as the pseudo-Mackay type clusters and the 
Bergman type clusters. It was observed that adjacent 
clusters were connected along either a two-fold or a three- 
fold symmetry axis of the reference icosahedron; a two- 
fold linkage connects clusters of the same kind and a 
three-fold linkage connects ones of different kinds. Af- 
ter the central atoms of all the clusters were identified, 
the missing atoms within each cluster were located by 
inspecting the charge density plot. Finally the refine- 
ment converged with the reliability index of R = 11.96 
(%) or R w = 12.13 (%). A detailed presentation of the 
refinement as well as a table of atomic parameters is pro- 
vided in Appendix <JB] The basic crystallographic data 
are summarized in Table [I] 



FIG. 5. The idealized templates for the two kinds of cluster, 
(a) M-cluster consists of a central site (Mo), a dodecahedral 
shell (M3) and a composite shell with an icosahedral subshell 
(M5) and an icosidodecahedral subshell (M2). (b) B-cluster 
consists of a central site (Bo), an icosahedral shell (B5) and a 
dodecahedral shell (B3). 




A. Cluster templates 

The geometrical templates for M- and B-clusters are 
shown in Figure [5] both have an atomic site at the cen- 
ter and two shells having the full icosahedral symmetry, 
532/m (ih)- The central sites are symbolized as Mo and 
Bo for M-cluster and B-cluster, respectively. The inner 
shell of M-cluster consists of twenty sites forming the ver- 
tices of a regular dodecahedron. This shell is symbolized 
as M3 because the relevant sites are on the three-fold 
rotation axes of the icosahedral point group; the same 
principle also applies when symbolizing the remaining 
shells. The outer shell of M-cluster is a composite of 
two subshells, twelve sites forming a regular icosahedron 
(M5) and thirty sites forming a regular icosidodecahedron 
(M2) . On the other hand, the inner shell of B-cluster con- 
sists of twelve sites forming a regular icosahedron (B5), 
while the outer shell twenty sites forming a regular do- 
decahedron (B3). 

Let us now define the six icosahedral basis vectors a.j 
(j = 1, 2, ... and 6) as 



ai a 2 a 3 a 4 a 5 a 6 

t 1 1 -r 
1 

-T 




(1) 



IV. CRYSTAL STRUCTURE 

The refined structure contains 4728 atomic sites per 
unit cell, wherein 204 sites are symmetrically unique. Im- 
portantly, the structure can be pictured as one formed 
through dense packing of clusters, which are allowed to 
overlap with each other across their peripheries. There 
is no need for glue atoms to fill in the gaps between 
the clusters, meaning each of the atomic sites belongs 
to at least one cluster. Two kinds of cluster called the 
pseudo-Mackay type and the Bergman type clusters are 
identified; these are henceforth referred to as M- and B- 
clusters, respectively. 



where r = (1 + \/5)/2 is the golden mean and I repre- 
sents an appropriate scale. The norm of the basis vec- 
tors is hereafter denoted a(:= | | ) . Note that the num- 
bering, j — 1, 2, ... and 6, of the six basis vectors is 
taken such that they are arranged concentrically around 
a three-fold rotation axis^. Taking the central site of 
a cluster template as the origin, the position vector x 
of every site within the cluster template is written as 
x = J2j c j a j =: [ c i c 2C3C4,C5Ce], in which the indices Cj 
(j = 1, 2, ... and 6) are all integers or half-integers. To 
be more specific, each of the shells, M^ and B K (x = 0, 2, 3 
or 5), is obtained as an orbit of its representative member 
with respect to the point group 532/m (Ih)', the repre- 
sentative indices and the radius are listed in Table [II] 



6 



TABLE II. The indices given for each shell represent the 
vector from the center to a representative site of the shell. 
Representative vectors for b- and c-linkages as well as the 
lattice translation vectors are also given for reference. The 
norm of each vector is given in the right column. Note that 
b — 2(r 3 /\/5) 1 ^ 2 a and b = y/3/2c. Negative indices are rep- 
resented by integers with top bars. 



snell 


indices 


norm 


Mo 


[000000] 





M 2 


[110000] 


6/t 


M 3 


[llllll]/2 


c/t 2 


M, 


[llllll]/2 


ra 


Bo 


[000000] 





B 3 


[lIllllj/2 


c/r 


Br, 


[100000] 


a (= 1.902£) 


b- linkage 


[110011] 


b (= 5.236£) 


c-linkage 


[111000] 


c (= 4.535^) 


Ri 


[604460] 


2fer 2 


R2 


[460046] 


2br 2 


R. 


[046604] 


26r 2 



TABLE III. Centers of the clusters within the unit cell ex- 
tracted from the refined parameters (Appendix fjBj) . 



site symbol atom Wyckoff indices node type 



B (l) 


Pdl 


8c 


[111000] 


(67)333 


d(J) 

B o 


Pd6 


8c 


[444111] 


(67)333 


B (3) 


Pdll 


24d 


[332001] 


(66)4322 


B (4) 


Pd24 


24d 


[443002] 


(67)333 


0(6) 

a o 


Pd37 


24d 


[312120] 


(76)433 


0(6) 

B o 


Pd50 


24d 


[423121] 


(76)433 


B<*> 


Pd63 


24d 


[533112] 


(67)333 


M« 


Cr/A176 


4a 


[000000] 


(68)0 


M< 2) 


Cr/A179 


4b 


[302230] 


(68)0 


M< 3) 


Fe/Pd82 


24d 


[432102] 


(57)3322 


M< 4) 


Fe/Pd90 


24d 


[412221] 


(57)332 


M< 5) 


Cr/A198 


24d 


[221001] 


(66)432 


M< 6) 


Fe/Pdl05 


24d 


[451024] 


(67)333 


M 7) 


Fell3 


24d 


[433011] 


(67)333 



B. Packing geometry 

In the refined structure (Appendix CBj). the unit cell 
accommodates a total of 264 clusters, which divide into 
128 M- and 136 B-clusters. In Table Hill the independent 
positions of the cluster centers are summarized. The cen- 
ters of the M-clusters (M ) are occupied mainly by Fe or 
Cr, although better fit is attained by occupying some 
of them partially with Pd or Al. The centers of the B- 
clusters (B ) are purely occupied by Pd. 

Adjacent clusters are mutually connected either along 



a two-fold or a three-fold axis with a distance of about 
7.7A or 6.7A, respectively; these are the shortest two 
distances between cluster centers. An analogous feature, 
with the distances being appropriately scaled, has been 
extensively discussed in the case of P-type icosahedral 
quasicrystals and their approximant o 8 ' 14 ^ 15 , in which only 
a single kind of cluster would come into play. The two- 
and three-fold linkages are called b- and c-linkages, re- 
spectively, whereas a packing of clusters with this prop- 
erty is called a bc-packing. The b- and c-linkages are in- 
dexed with integers and their representatives are given in 
Table [Hi the respective norms are denoted b and c. The 
skeleton of the crystal structure can thus be described 
as a network of nodes, which correspond to the cluster 
centers, connected through b- and c-linkages. 

F-type ordering (or F-centering) in icosahedral qua- 
sicrystals as well as their approximants occurs as the even 
and odd parities of the nodes in the relevant bc-packing 
are differentiated. Here, the parity of each node is de- 
fined as that of the sum of the relevant indices. Note in 
TablellTlthat each b- linkage connects a pair of nodes hav- 
ing the same parity, while that each c-linkage connects a 
pair of nodes having different parities. Therefore, if the 
two kinds of cluster, M and B, are the entities that dif- 
ferentiate the two subsets of the nodes of a bc-packing, 
it follows that every b-linkage connects the same kind of 
cluster (M-M or B-B), while that every c-linkage connects 
different kinds of cluster (M-B). The three combinations 
for an adjacent pair of clusters are depicted in Figure [SJ 

The coordinates of the cluster centers given in the re- 
fined atomic parameters (Appendix fjB| can be used to 
enumerate which of the b- and c-linkages connect every 
site to its adjacent neighbors. The task of obtaining the 
indices for every cluster center is straightforward. First 
we set the indices of the atomic site Cr/A176 lying at 
the origin (0,0,0), which is the center of an M-cluster, 
to [000000]. Then we recursively trace linkages to ob- 
tain the indices of adjacent cluster centers until all the 
cluster centers within the unit cell as well as the lattice 
translation vectors are identified. The resulting lattice 
translation vectors Rj (j = 1, 2 and 3) are presented in 
Table |TTJ Since the lattice constant, IRjl = 2br 2 should 
equal the experimental value of a/ at = 40.54A, the basic 
parameter b is immediately evaluated as 7. 74 A. 

Now we are equipped with an appropriate scale for ex- 
amining the inter-atomic distances. Let us first note that 
the ideal edge length 6/r 3 of the inner dodecahedral shell 
of M-cluster (M 3 ; FigureEJa)) is evaluated to be 1.83A, 
and that it is unrealistically short for an interatomic dis- 
tance. A natural consequence of this is that the relevant 
shells cannot be occupied by more than 8 atoms at the 
same time. Next, note that the outer shell of M-cluster 
(including the two subshells) has two kinds of edges. The 
edges connecting the adjacent pairs in the M 2 subshell 
are parallel to the two-fold axes and have a length of 
b' := b/r 2 = 2.958A, while those connecting the M2 sub- 
shell with the M5 subshell are parallel to the three-fold 
axes and have a length of c' := c/t 2 = 2.562A. The latter 
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FIG. 6. The three different combinations for an adjacent pair 
of clusters: M-M (left), B-B (middle) and M-B (right). Two 
clusters of the same kind are connected via b-linkage, while 
those of the two different kinds are connected via c-linkage. 
The sites in the constituent cluster templates are depicted 
as spheres. The drawing plane is perpendicular to a 2-fold 
symmetry axis. 




two edge lengths are reasonable for interatomic distances. 
Similarly, the edges connecting the adjacent pairs in the 
inner (B5) or the outer (B3) shell of B-cluster are parallel 
to the two- fold axes and have a length of b' . Interest- 
ingly, the two kinds of linkages, namely b' and c', form 
the majority of the interatomic linkages in the idealized 
construction. This applies not only to the closest inter- 
atomic linkages within each shell of a cluster but also to 
those connecting between the inner and the outer shells. 
It follows that the closest distances are b' for the pairs, 
M3-M3, M 2 -M 2 , M3-M5, B5-B5 and B3-B3, while c' for 
the pairs M2-M5, M3-M2 and B5-B3. Here we have in- 
cluded the pair M3-M3 because the closest interatomic 
distance that is allowed within the partially occupied in- 
ner shell (M3) of an M-cluster is b'. 

Two M-clusters connected through b-linkage interpen- 
etrate into each other (Figure^ left) with the overlap be- 
tween the outer polyhedra being a flat hexagonal bipyra- 
mid. The two tips of the bipyramid belong to the M2 
shells and are very close to each other (1.828A), so that 
they cannot be occupied at the same time. In the re- 
finement, such positions have been fitted as splitting po- 
sitions for a single Al atom. The intersection between 
the two M-clusters consists of the six sites forming the 
base hexagon of the hexagonal bipyramid; two of them 
belong to the M5 subshells of the two M-clusters, while 
four of them to the M2 subshells. On the other hand, two 
B-clusters connected through b-linkage share an edge be- 
tween their outer shells (B 3 , Figure [6l middle). In both 
cases, the intersection between the two clusters does not 
involve the inner shells. 

The situation is somewhat more intricate in the case 
of an M-B pair connected through c-linkage. Observe in 
Figure E] (right) that the outer shell (B3) of the B-cluster 
penetrates to the inner shell (M 3 ) of the M-cluster, where 
a site is shared by these shells. Similarly, the outer sub- 
shell (M2) of the M-cluster penetrates to the inner shell 
(B5) of the B-cluster, where a triangular face is shared 
by these shells. One also finds that these clusters inter- 



FIG. 7. The packing of clusters. The square drawn with white 
dashed lines indicate a face of the cubic unit cell, whereas 
the horizontal and vertical edges correspond to the primitive 
lattice vectors, Ri and R2, respectively. 




sects with each other in their outer (sub) shells M5 and 
B 3 , too, where there are three common sites between the 
two. 

The idealized atomic positions can be obtained by 
replicating the relevant cluster templates at the posi- 
tions given in Table Mil The arrangement of clusters in 
a unit cell is illustrated in Figure [3 As we have just 
seen, an atomic position can be shared by two or more 
cluster shells that intersect with each other. Therefore, 
each individual site can be characterized by the set of 
the cluster shells to which it belongs. Take, for instance, 
a site in the inner shell (M3) of an M-cluster belong- 
ing also to the outer shell (B3) of an adjacent B-cluster. 
We simply assign to this site a class, (M3,B3). More 
generally, if an atomic site belongs to n cluster shells, 
X , X 2 , X n , then the class associated with this site 
would be (A 1 , A 2 , A"). Bear in mind that only the 
combination of the cluster shells matters here, so that 
the order of the shell symbols in the angle brackets are 
irrelevant. It turns out that the idealized structure of 
the present crystal contains in total 16 different classes 
of atomic sites, which are listed in Table HVl 

It is obvious from the above argument that the central 
site of each cluster cannot be shared by any other cluster. 
Hence, the symbols for the cluster centers are (Mo) and 
(Bo). In addition, the crystal structure includes three 
more classes of unshared sites, which are (M.2),(Ms) and 
(B 5 ). Importantly, the former two bear distinct roles 
in the physical construction of the structure: (i) (M2), 
the splitting positions given at the tip of the overlap 
hexagonal bipyramid associated with each interpenetrat- 
ing pair of M-clusters and (ii) (M3), the vacant sites in 
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TABLE IV. The 16 different classes of idealized atomic sites. 
Their number frequencies per unit cell, presented in the sec- 
ond column, sum up to a total of 4680. The average composi- 
tion (in percentage) of each class is given in the last column, 
where 'Vc' stands for vacancy. 



site class 


number 


average composition 


(M ) 


128 (= F) 


Al 6 Pd 4 (Cr,Fe) 9 o 


/TV IT \ 

<M 2 ) 


720 (= G) 


Al 5 oVc 50 


(M 3 ) 


1680 (= H) 


Vcioo 


<Bo) 


136 (= I) 


Pdioo 


/T> \ 

( B 5> 


24 (= J) 


AI100 


(M 2 ,M 2 ) 


72 (= K) 


AI100 


<M 2 ,B 5 ) 


576 (= L) 


AI100 


(M 3 ,B 3 ) 


880 (= M) 


AlgiPdi(Cr,Fe)i 7 Vc : 


<M 5 ,B 3 } 


24 (= N) 


Pdioo 


<B 3 ,B 3 ) 


72 (= O) 


Pdioo 


(M 2 ,M 2 ,B 5 ) 


1032 (= P) 


AI100 


(M 2 ,M 2 ,M 2 ) 


112 (= Q) 


AI100 


(M 5 , M 5 , M 5 ) 


24 (= R) 


Ali(Cr,Fe) 9 9 


(M 5 ,B 3 ,B 3 ) 


408 (= S) 


Pd 98 (Cr,Fe) 2 


(M 5 ,M 5 ,B 3 ,B 3 } 


384 (= T) 


Pd 6 5(Cr,Fe) 35 


<M5,M 5 ,M 5 ,B 3 } 


88 (= U) 


Pd 3 (Cr,Fe) 97 



the M3 shells of the M-clusters. Remember that the inner 
shell of M-cluster (M3) cannot accommodate more than 8 
atoms due to the constraint imposed by the short nearest- 
neighbor distance. And it is clearly demonstrated by our 
structure analysis that the only occupied M3 sites are 
those represented by the symbol (M 3 ,B 3 ). The last un- 
shared class (B5) does not seem to differ from the other 
types of B5 sites, as all of them are fully occupied by Al. 

Correlations are further found between the local com- 
positions at individual atomic sites and the site classes. 
Take, for instance, the B3 shell, which involves six dif- 
ferent site classes (Table IIV[) . Whereas the sites sym- 
bolized as (M3,B3) are primarily occupied by Al, the 
sites in the remaining five classes, (M 5 ,B 3 ), (B 3 ,B 3 ), 
(M 5 ,B 3 ,B 3 ), (M 5 ,M 5 ,B3,B 3 ) and (M 5 ,M 5 ,M 5 ,B 3 ), are oc- 
cupied by markedly heavier elements in different degrees 
depending on their site classes. Similar observations sug- 
gest that the physical nature of an individual atomic site 
is profoundly affected by the manner how clusters inter- 
sect there. Therefore, the two kinds of clusters not only 
provide a handy means of describing the complex struc- 
ture but also serve as the true physical units which play 
a significant role in the formation of the structure. This 
underlying basic idea may also apply to a more general 
class of Al-based alloys which exhibit F-type icosahedral 
ordering. 

The packing of the cluster templates as described 
above generates in total 4680 idealized atomic positions 
(excluding the vacant (M3) sites), which agree remark- 
ably well with the refined atomic positions, although in 
the refinement we included 48 additional sites which are 
needed to explain minor irregularities in the real mate- 
rial. Importantly, the gaps between the clusters are small 



enough and do not accommodate any additional atoms 
(called glue atoms). In Subsection ^VCi this feature is 
described from a somewhat different viewpoint. 



C. Canonical cell tiling with F-type ordering 

Recall that the atomic arrangement and the chemical 
compositions within each individual cluster are subject 
to constraints posed by the existence of adjacent clus- 
ters (Subsection ^IVBI) . These geometrical constraints 
may hardly allow other combinations of clusters than the 
three cases shown in Figure |6l thus enforcing a perfect F- 
type ordering. The F-type icosahedral ordering observed 
in more general Al-based alloys could also have a sim- 
ilar geometrical origin, where an uneven distribution of 
atomic species on different sublattices is understood to 
be a secondary consequence. 

The skeletal structure of the present approximant is 
described as the bc-packing, in which the cluster centers 
are represented as the nodes (Table IIIII) . The nodes can 
be connected to each other through edges of the two kinds 
(b- and c-linkages; Figure HJa)), whereas the edges form 
three kinds of polygons: an isosceles triangle formed by 
a b-linkage and two c-linkages (X-face; point symmetry, 
m), an equilateral triangle formed by three b-linkages 
(Y-face; 3m) and a rectangle formed by two b-linkages 
and two c-linkages (Z-face; 2/m); Figure [SJ^b) . These 
polygons are further found to be the faces of four kinds of 
polyhedra called the canonical cellA a tetrahedron with 
four X- faces (A-cell; 2m), a pyramid with three X- faces, 
one Y-face and one Z-face (B-cell; m) , a tetrahedron with 
three X-faces and one Y-face (C-cell; 3m) and a trigonal 
prism with two Y- faces and three Z-faces (D-cell; 3m); 
Figure [H](c). The present bc-packing thus proves to be 
represented as a periodic tiling of space called a canonical 
cell tiling (CCT); that is, the whole space is divided into 
pieces congruent to the canonical cells (Figure [9]) . 

So far, suitable atomic decorations of periodic CCT's 
have been commonly used to describe a variety of approx- 
imant phases to P-type icosahedral quasicrystal o 8 i 15 . On 
the other hand, the present compound demonstrates for 
the first time that the CCT construction extends natu- 
rally to approximant phases to F-type icosahedral qua- 
sicrystals. Here, the two parities of the vertices need to 
be clearly distinguished to account for the F-type order- 
ing; this is done by representing even (resp. odd) vertices 
with white (resp. black) spheres in Figure |H1 The geo- 
metrical composition for the present CCT in particular is 
summarized in Table |Vj There are 1 1 congruence classes 
of objects, some of which are further divided into two 
subclasses (I and II) owing to the parities of the nodes, 
leading to the 18 object classes in total. Remember that 
X-face as well as A-, B-, C- or D-cell has more than one 
unique corners (see the corner symbols in Table |V|) . for 
which the parities need to be specified explicitly. It fol- 
lows that there are two ways in coloring the vertices of 
each kind of cell, leading to two different configurations 
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FIG. 8. The canonical cells and their geometrical components, 
(a) The two kinds of linkages, b and c, represented as thick 
and double bars, respectively, (b) The three kinds of faces, 
X, Y and Z, whose sides are distinguishable in terms of the 
icosahedral symmetry; where '+' sign is on one side, ' — ' sign is 
on the other side. The reverse sides of the faces are shown by 
rotating them by 180 degrees around the vertical dashed line, 
(c) The four canonical cells, A, B, C and D, and a trigonal 
anti-prism unit, B2, which can be divided into two B-cells in 
three ways. The unique corners are indicated by the corner 
symbols, while the numbers on the edges correspond to the 
following dihedral angles: (1) 7r/2, (2) 7r/3, (3) ir — e, (4) 7r — 77, 
(5) 17, (6) 7T- 277, (7) (7r-e)/2, (8) e, (9) 77, (10) tt/2 and (11) 
tt/3. 




of clusters and thus to two different atomic arrangements 
within the cell (in the sense that the atomic species are 
disregarded). 

The present argument extends naturally to arbitrary 
CCT's, whereby a number of hypothetical approximants 
with F-type ordering can be constructed. From this 
broader perspective, it is worthwhile making a general 
consideration on the statistics of the geometrical objects 
as well as that of the atomic sites. Still, readers who 
might want to avoid handling the mathematical deriva- 
tions may safely skip the rest of this subsection and go 
to Section JV] 

In the following argument, the number density of an 
arbitrary object type O among the 18 basic object types 
is denoted n(0). In particular, the number densities of 



FIG. 9. (color) A top view of a cleaved surface of the present 
CCT parallel to the (001) plane; the foremost cells are cho- 
sen somewhat arbitrarily. The spheres representing the even 
(resp. odd) nodes are colored white (resp. black). The white 
dashed lines indicate a square face of the cubic unit cell, while 
the horizontal and vertical edges correspond to the primitive 
lattice vectors, Ri and R2, respectively. At the lower left cor- 
ner of the square lies an even node corresponding to the origin. 
The four kinds of canonical cells are colored yellow (A-cell), 
red (B-cell), blue (C-cell) and green (D-cell), respectively. 




TABLE V. The number frequencies of the geometrical objects 
contained in a single unit cell of the present CCT. The 11 con- 
gruent types can be further divided into subclasses according 
to the parities of the nodes as specified in the parentheses us- 
ing the symbols, + (for even) and — (for odd), leading to 18 
basic object types. Note that for the faces (X - Z) as well as 
the cells (A - D), the parity symbols are arranged in the same 
order as the corner symbols referenced in the second column. 
The object type B2 can be divided into two B-cells in three 
different ways, and we do not take it as a basic object type. 



object 


type 








(I) 




(II) 


node 


128 (+) 




136 (-) 


edge b 


360 (+ + 


) 


432 (- -) 


edge c 


* 880 (+ 


-) 




face X (X X a Xi,) 


1032 (+ - 


«— ) 


1176 (- - +) 


face Y (Y a Y a Y a ) 


112 (+ + 


+) 


176 ( ) 


face Z (Z a Z a Z a Z a ) 


* 176 (+ 


+ --) 




cell A (AaAaAfcAi,) 


384 (+ + 


--) 


336 (-- + +) 


cell B (BaBaBtBfcBc) 


104 (- - 


+ + +) 


152 (+ + 


cell B 2 (B C B C B C B C B C B C ) 


* 80 (+ -I 


- H 


-) 


cell C (CaCsCfiCb) 


88 (- + + +) 


168 (+ ) 


cell D (DaDaDaDjDfcDj,) 


8 (+ + + 


) 


24 ( + + 
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the eight types of cells are of fundamental importance, so 
that we introduce the following parameters: a\ — n(Aj), 
a u = n(A n ), (3i = n(Bi), /3 n = n(B n ), 71 = ra(Ci), 
7n = n(Cu), Si — n(Di) and Su = n(Dn). The following 
sums are also used whenever it is convenient to do so: 



a:=a>i+aii, /3 := j3\ + /3 n , 
7 : = 7i +7ii) <5 := Si + 8 U . 



(2) 



By definition, the sum rule for the cell volumes is written 



VAa + vb/3 + vcl + vdS = 1, 



(3) 



where v A = b 3 /12, v B = Vbb 3 /12, v c = V$b 3 /2A and 
vo = 3b 3 /8 are the volumes of A-, B-, C- and D-cells, 
respectively. In addition to this, these statistical param- 
eters hold a set of universal equations due to geometrical 
constraints imposed by the shapes of the cells. 

Let us first consider the matching constraints across 
the faces. Take, for instance, a cell that has an Xi face 
with its + side facing outward, then another cell with an 
Xi face with its — side facing outward can only match the 
former cell. Since the same requirement is fulfilled across 
any of the Xi faces in the structure, the + and — sides of 
the faces of this type must appear equal number of times. 
By counting the two distinct sides of Xi faces for each 
type of cells, Eq.Q can be readily verified. Furthermore, 
similar considerations when applied to the Xn, Yi and 
Y11 faces lead to Eqs.©, ((SJ) and ([7]), respectively. 



2ai + 3 7 i = 2a„ + 2/3i + Ai, 
2aii + 3 7 „ = 2ai + 2/3„ + A, 
/3i + 61 = 71 + Sn, 
fin + $11 = 711 + Si. 



(4) 
(5) 
(6) 
(7) 



Now we show that the matching constraints around 
the edges will lead to two more equations. Note that the 
dihedral angles contained in a B- or C-cell (see Figure 
IHfc)) involve two irrational constant, e and 77, defined 
by 8 



e = arccos(i) w 4tt(0.1049), 



77 ee arccos (—=) « 4tt(0.0915). 
V6 



(8) 
(9) 



Here each cell contributes a fixed amount of dihedral an- 
gles associated with a fixed edge type. Hence, the con- 
tributions from all the cells sum up to the total number 
of the relevant edges times 2tt; that is, 



2nn(bi 



■vPii 



{w/2)a + (3tt - 4r?)/3i 
+37771 + (37r/2)5, 
2vm(b„) = (7r/2)a + r)fr + (3tt - 4t7)/3„ 
+377711 + (37r/2)<5, 
27rn(c) = (47r/3)a + 3(tt - e)/3 + 3e7 + ttS. 



(10) 

(11) 
(12) 



Observe that in Eqs. (fl"U)) - (|12[) terms containing e and 77 
do not cancel out spontaneously. Hence, in order to avoid 



unphysical consequences that is the number of edges of 
each type would be irrational, we need to constrain the 
coefficients for e and 77 to be zero. Therefore, the follow- 
ing equations need to be satisfied; 



+ 371 

-3/3 + 3 7 



A - 4/3 n + 3 7n - 0, 



(13) 



One can readily check that Eqs.Q-© and (fTB")) can be 
reduced to the following four universal equations, 



where 



= 7, Aa = -Afl 
3A 7 = 5A/3, 3AS = A/3, 



Aa := ai - a n , A/3 := /3/ - j3n, 
A7 := 7/ - 7//, AS := (5/ - 5 H . 



(14) 



(15) 



It follows from Eqs.© and (fl4l that there are only 
three independent degrees of freedom to determine the 
statistical properties of any CCT with F-type ordering. 
Following ref|8|, an independent parameter \x is defined 
as the volume fraction occupied by A- and D-cells: 



/j := VAa + v D S, 



(16) 



while another £ is defined so that gives the volume 
fraction occupied by D-cells: 



C :=v d S/ii. 



(17) 



It is readily shown that these two parameters determine 
the number densities of the canonical cells through 8 



a = 12/41 - C)/6 3 , 

/3EE 7 = (8/V5)(l-/i)/& 3 , 

S = (8/3R/& 3 - 



(18) 
(19) 
(20) 



The last independent parameter v, which is necessary if 
F-type ordering is considered, can be defined as 



v := A/3//3. 



(21) 



As soon as the statistics is fixed for the cells, it is a 
straightforward task to enumerate the objects of lower 
dimensionalities (faces, edges and nodes) by taking the 
contributions from all the cells. The edges have already 
been enumerated via Eqs.([TU l) - p^)) . while for enumerat- 
ing the faces the double counting of each face must be 
taken care of. The nodes are enumerated by taking the 
sum of the solid angles associated with the relevant cor- 
ners of the cells^ and by dividing the results by Air. The 
resulting formulae are 

77(nodci) = a/12 + /3/4 + A/3/12 + S/A, (22) 

n(node„) = a/12 + /3/4 - A/3/12 + S/A, (23) 

n(bi) = a/4 + 3(/3 + A/3)/ A + 36/4, (24) 

n(b„) = a/4 + 3(/3 - A/3)/4 + 36/4, (25) 

n(c) = 2a/3 + 3/3/2 + S/2, (26) 

n(Xi) = a + 3(/3 + A/3)/2, (27) 
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n(Xn) = a + 308 - A/3)/2, (28) 

n(Yi) = /3/2 + 2A i 8/3 + 5/2, (29) 

n(Y„) =,3/2-2A/3/3 + 5/2, (30) 

n(Z) = (/? + 3<5)/2. (31) 



The numbers of objects listed in Table IVl are the frequen- 
cies within a single unit cell. These values can be read- 
ily compared with the above formulae by re-interpreting 
the variables as the frequencies per unit cell, that is, 
aQ = 720, (3ft = 256, AfiQ = -48 and 5ft = 32. 

Each of the above 18 object types offers a possible 
local atomic configuration through replicating the rel- 
evant cluster templates. Therefore, the above statis- 
tics uniquely determines the number densities / = F/fl, 
g = G/Cl, u — U/Sl of the 16 types of atomic sites; 
see Table ITVl One can readily prove the following general 
formulae (see Appendix [JXJ: 



/ = 


l/12a + 1/4/3 + 1/12A/3 + 1/4(5, 


(32) 


9 = 


1/2q + 3/2/3 + 3/2A/3 + 3/25, 


(33) 


h = 


a + 7/2/3 + 5/3A/3 + 9/25, 


(34) 


i = 


l/12a + l/4p- 1/12 A/3 + 1/45, 


(35) 


i = 


1/2AP + 3/2S, 


(36) 


k = 


-1/2A/3 + 3/25, 


(37) 


I = 


3/2/3 -3A/3 + 3/25, 


(38) 


m = 


2/3a + 3/2/3+ 1/25, 


(39) 


n = 


l/2A^ + 3/25, 


(40) 


o = 


-1/2A/3 + 3/25, 


(41) 


P = 


a + 3/2/3 + 3/2A/3, 


(42) 


q = 


1/2/3 + 2/3A/3+ 1/25, 


(43) 


r = 


-1/6A/3+ 1/25, 


(44) 


s = 


3/2/3- 1/2A/3, 


(45) 


t = 


l/2a- 1/2A/3, 


(46) 


u = 


1/2/3 + 5/6A/3. 


(47) 



Remembering that each individual site exhibits a chemi- 
cal composition that is substantially correlated with the 
site class, the derived statistics could be useful when esti- 
mating the chemical composition of a hypothetical crys- 
tal that is derived from a CCT. This may facilitate an 
experimental search for unknown approximants in closely 
related alloy systems. 

V. DISCUSSION 

The new model presented in Section ^IVI allows a con- 
cise description of the extremely complex material, offer- 
ing the possibility that the stabilization mechanisms as 
well as the physical properties will be clarified on a struc- 
ture basis. In particular, the model will lead us to revised 
views on bondings, defects and dynamics in related mate- 
rials. At present, our model cannot be extended naively 
to icosahedral quasicrystals because there is currently no 
proof as to whether a CCT— can be made quasiperiodic 



nor having the icosahedral symmetry 16 . One can never- 
theless gain an important insight into the quasicrystallinc 
state by evaluating how close (or how far) the given crys- 
tal lies to (or from) an idealized quasicrystal. The fol- 
lowing subsections are devoted to discussions on some of 
these topics. 

A. Superlattice ordering 

In the standard six-dimensional formalism for icosahe- 
dral quasicrystals, the three-dimensional physical space 
E is thought of as a subspace of a six-dimensional hyper- 
space E, whereas the orthogonal complement to the phys- 
ical space is called the perpendicular space Ej_. The 
hyper-space is the direct sum of the two subspaces; 
i.e., E = E © Ej_. Accordingly the basis vectors a.j 
(e E) defined in Eq.(fT]) are lifted into six dimensions 
via &j = ca.j), where c(s K) is an arbitrary non-zero 
factor and aj- (e Ej_) are defined by 

( a l a 2 a 3 a 4 a 5 a S ) 

/ 1 -r -r -1 \ 
:=£ -r 1 -r -1 . (48) 

y o — t i — i o — t J 

If c is taken to be unity, the six-dimensional basis vectors 
&j (e E) form an orthogonal basis set that generates a 
six-dimensional hyper-cubic lattice, Lp, where the sub- 
script P stands for the primitive lattice (P-type). The 
norm of the six-dimensional basis vectors are hereafter 
denoted a(:= |eLj| = \/2a), which is evaluated to be 3.977 
A from the present approximant. 

Provided that we have a structure in which the ar- 
rangement of M- and B-clusters is subject to the local 
rules presented in Subsection ^IVBl each cluster is cen- 
tered at a position that is represented as 

x = mai + n 2 a 2 + n 3 a 3 + n 4 a 4 + n 5 a 5 + n 6 a 6 
= [nin 2 n 3 ri4n 5 n 6 }, (49) 

where rij (J = 1, 2, ... and 6) are the indices, i.e., the 
integer coefficients that are determined uniquely if the 
origin is taken at one of the cluster centers. The kind of 
the cluster (resp. M or B) depends on the parity of x 
(resp. even or odd), where the parity is defined as that 
of the sum Y^, , — n j ■ 

The three-dimensional vector x can be lifted to a six- 
dimensional counterpart x € Lp simply by replacing &j 
in Eq. (|4"9"]) with a,-. Then, x is the image of x through the 
orthogonal projection of E onto E. The image x- 1 of x in 
Ej^ is definable in a similar way. It is readily understood 
that the three-dimensional vectors x and x^ maintain 
a one-to-one correspondence (bijection) with each other, 
where x 1 - is called the conjugate image of x. The parities 
of x, x and x are equal. The even subset (consisting of 
the vertices of even parity) in Lp forms a sublattice with 
an index of 2 and it gives an F-type (i.e., F-centered) 



12 



hyper-cubic lattice Lp, having a lattice constant of of := 
2a. 

The perfect icosahedral symmetry of an ideal qua- 
sicrystal requires that the conjugate images of all the 
cluster centers be bounded within a finite domain W 
(c Ex) called a window, which is also called an atomic 
surface or an acceptance domain. If this is the case, 
the lifted coordinates x are distributed along the three- 
dimensional cut space E + w in E, where w (G Ej_) is 
the center of mass of the window. This can be put also 
as x 1 - wo, which is a loose but useful expression in a 
context where the deviation from the equality is unim- 
portant. 

If the cut space is taken to be inclined linearly against 
E, the coordinates of the cluster centers would behave as 
x 1 - ~ A • x + wo , in which A is a 3 x 3 matrix called the 
linear phason tensor. An approximant in particular has 
a cut space that is parallel to one of the lattice planes 
of Lp , such that any lattice translation vector R for the 
approximant satisfies R^ = A • R, where R^ is the con- 
jugate image of R. Therefore, the task of evaluating A 
for the present approximant can resort to the primitive 
lattice translation vectors given in Table HH reading 

A = ( R^ R^ - R^j" ) ( Ri R2 R3 ) 

29-13V5 T 7 , s 

= 2"-^ I = -r- 7 I, (50) 

where I is the 3x3 identity matrix. The smallness of 
the linear phason strain implies that the crystal is rather 
close to an ideal quasicrystal. 

The present approximant can indeed be identified as 
a 2 x 2 x 2 superstructure of the cubic 3/2 approximant 
in the series of rational approximantions of icosahedral 
quasicrystalinity. The ordinary cubic 3/2 approximant is 
constructed from a CCT called the 3/2 packing 8 -, which 
has a lattice constant of br 2 ; then, a hypothetical crys- 
tal structure illustrated in Figure [TU1 is obtained using the 
decoration rules of Section SjlVl Observe a marked resem- 
blance of this unit cell to each 1/8 block (half cube) cut 
out from the cubic unit cell of the superstructure (Figure 
[7]) . The primitive lattice translation vectors of the hypo- 
thetical 3/2 approximant are R^- = R^/2 with j = 1, 2 
and 3. Note that the linear phason tensor is unchanged 
from that given in Eq. (1501) because the corresponding in- 
dices (even numbers) for Rj are simply halved for R^ . 
The space group Pa3 of the 3/2 packing 8 is reduced 
to P2i3 because the atomic decoration distinguishes the 
even and odd parities of the nodes. Accordingly, the hy- 
pothetical 3/2 approximant should have the enantiomor- 
phic (or chiral) point group 23. When the superlattice 
ordering is induced, the space group Pa3> including the 
center of symmetry is restored. 

It turns out that not only rearrangements of clusters 
but also additions of extra clusters close to the bound- 
aries between half cubes are necessary in realizing the 
superstructure. The latter can be checked when the pack- 
ing densities of the two structures are compared. Observe 



FIG. 10. The packing of clusters for a hypothetical cubic 
3/2 approximant. The square drawn with white dashed lines 
indicates a face of the cubic unit cell, whereas the horizontal 
and vertical edges correspond to the primitive lattice vectors, 
R' x = Ri/2 and R2 = R2/2, respectively. 




that a single unit cell of the 3/2 approximant with a vol- 
ume of b 3 r 6 (=: ^3/2) contains 16 M- and 16 B-clusters, 
while that one of the present superstructure with a vol- 
ume of = 8^3/2 contains 128 M- and 136 B-clusters. 
Eight additional B-clusters (i.e., 8 = 136 — 16 x 8) are 
included in the unit cell, causing an increase in the num- 
ber density of clusters from 1.78/fo 3 to 1.84/6 3 . The 
present superstructure gives a relatively dense packing 
as compared to most of the simplest CCT's considered 
by Henley 8 -; the only two with higher densities are the 
cubic 1/1 and 2/1 packings with the number densities of 
clusters being 2/b 3 and 1.89/& 3 , respectively. 

A clear contrast between the basic skeleton of the 
present approximant (superstructure) and that of the 
3/2 approximant can be demonstrated with their modi- 
fied conjugate images. Generally speaking, the conjugate 
images of the cluster centers would be unbounded if there 
exists a non-zero linear phason strain, i.e., A ^ O (O, 
the 3x3 zero matrix), because x in x x ~ A • x + wo 
runs across the physical space E. Hence, one needs to 
consider the modified conjugate images 

X mod := X± — A • X, (51) 

which behave as x^ od ~ Wo. In Figure fTTl the modified 
conjugate images of the nodes of the two skeletons are de- 
picted as the white and black spheres, which correspond 
to the centers of M- and B-clusters, respectively. The 
truncated octahedron & depicted in each panel (centered 
at the center of mass of the modified conjugate images of 
the nodes) represents to the modified window associated 
with the 3/2 packing. All the vertices and the hexagonal 
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FIG. 11. The modified conjugate images of the cluster cen- 
ters for the ordinary 3/2 packing (left) and for the 2x2x2 
superstructure (right). The black (resp. white) spheres rep- 
resent the images of the centers of the B-clusters (resp. the 
M-clusters). The weight (or occupancy) of each modified 
conjuage image is represented by the volume of the relevant 
sphere. The truncated octahedron G shown in each panel in- 
dicates the outer boundary of the fundamental domain in Ex; 
see text for the definition. 
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face centers of 6 are obtained as the modified versions 
(via Eq.(|5ip) of the 32 vertices of a rhombic triaconta- 
hedron ^2 (c EjJ, which is the window for the 12-fold 
vertices in the three-dimensional Penrose tiling; for the 
definition of , refer to the original article^ and Figure 
[9](a) therein. We find in Figure [TTT left) that the modi- 
fied conjugate images of the nodes in the 3/2 packing 
are well confined within G, where every point receives 
one node per unit cell. On the other hand, those for 
the 2x2x2 superstructure have wider spread than G 
(Figure [TTJ right). Importantly, the images lying close to 
the center of G receive full number of nodes (8 per unit 
cell), whereas those in the peripheral region receive less 
(<8). This indicates clearly that the windows associated 
with the relevant six-dimensional lattice Lp should be 
non-uniform. 

Relating the present superstructure in E (3D) to one 
in E (6D) is of significant importance in arguing fur- 
ther possibilities that superlattice ordering be observed 
in other approximants as well as quasicrystals. Indeed, 
it has been reported by Ishimasa£ that a superlattice 
ordering was exhibited by an icosahedral Al-Pd-Mn al- 
loy after being annealed at 700° C. Performing a careful 
analysis of diffraction patterns, Ishimasa concluded that 
the Bravais class for the superlattice was P-type, while 
the corresponding six-dimensional lattice constant was 
ap = 20.88lA= rap, where r is the golden mean and 
ftp = 12.901 A is the conventional six- dimensional lat- 
tice constant for the F-type icosahedral phase Al-Pd-Mn. 
One can check that the equation ap = rap holds, where 
dp is the lattice constant of Lp ; hence, the conventional 
lattice constant is r times that determined in the present 
study. The discrepancy is however non-essential since 
there exists an arbitrariness in choosing the basis vectors 
&j in E up to the multiplication of any integral power of r 
if only given the diffraction module, because of the r scal- 



ing invariance of the F-type icosahedral module and the 
corresponding I-type Fourier module^—. It is only be- 
cause the relevant scale for the cluster packing has now 
been fixed (Subsection § i]IVB[) that we prefer to use a 
more proper lattice constant dp instead of the conven- 
tional one ap. The lattice constant for the P-type super- 
lattice ordering in the quasicrystal is therefore written 
as ap = r 2 dp. Remember that the use of the updated 
lattice constant would affect the order of rational ap- 
proximation; e.g., a conventional cubic 2/1 approximant 
should now be called a cubic 3/2 approximant. 

The simplest P-type sublattice (or superlattice) of Lp 
is obtained just by taking the subset whose members have 

only even indices, i.e., Lp -* := 2Lp={ [711712713714715716]; 
nj = (mod 2) for j — 1, 2, ... and 6}. Then it is a sub- 
lattice of Lp with an index of 32 and its six-dimensional 
lattice constant is given by a p ^ — dp. Obviously, L p ^ 
fails to represent the superlattice ordering in the qua- 
sicrystal because a P °' ^ ap(— T 2 dp). In order to obtain 
the proper sublattice of Lp to describe the superlattice 
ordering of the quasicrystal, we need first to re-index Lp 
prior to taking the subset with only even indices. The re- 
indexing is done with respect to the re-scaled basis vec- 
tors aj 2 ' = (r 2 aj, c(— l/r) 2 a^), which can be obtained 
from the original basis vectors §Lj through 



5(2) s (2) ~(2) ~(2) ~(2) ~(2) 



3.A &c 



a 



(a^a! 
= (&i a 2 a 3 a 4 a 5 a 6 ) M 2 



(52) 



where six components for &j and aj 2 ' are all aligned 
column-wise and M represents the uni-modular r-scale 
transformation matrix 20 . 



M 



(53) 



/I 1 1 1 -1 -1 ^ 
1 1 1-11-1 

1 111-1-11 

2 1-1-11-1 -1 
-1 1-1-11 -1 

V -1 -1 1 -1-1 1 / 

(2) 

Note that the re-scaled basis vectors ; are not orthog- 
onal to each other if c = 1. Nonetheless, since the scale 
factor c can be arbitrarily chosen without affecting the 
projected images in E, one can choose c = r 4 so that the 

re-scaled basis vectors are orthogonal (i.e., a^p = T 2 &j) 

~ (2) 

and the sublattice L P with even indices would regain 
the hyper-cubic symmetry with a six-dimensional lattice 

constant of r 2 ap. This proves that Ishimasa's superlat- 

~ (2) 

tice ordering is associated with the Bravais lattice L p . 

There can be another kind of P-type sublattice L p 1 ' ) de- 
fined after re-indexing with respect to the re-scaled basis 
vectors a^- = (ra,-, c(— l/r)a+), which are obtained via 



sM sW s« ;(!) 5« s« 



(ai 



(aia 2 a3a4a 5 a 6 )M. 



(54) 
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When the perpendicular space Ej_ is appropriately scaled 
(i.e., c = — < ); the six-dimensional lattice constant for 

Lp** would be ra-p. 

The three sublattices, h p °\ hp and Lp-*, of Lp are 
distinct from each other in terms of translational sym- 
metry. One needs to clarify whether the superlattice 
ordering in the present approximant can be associated 
with the same sublattice L p ; as in the case of the qua- 
sicrystal. In order to simplify the argument, let us focus 
our attention on the M-cluster whose center lies at the 
body center B = [555111] of the conventional unit cell. 
Obviously, B is inequivalent to V = [000000], i.e., the 
vertex of the unit cell, through translation because the 
crystal is primitive cubic. The latter condition would 
be satisfied if the relevant superlattice ordering in six di- 
mensions were described by Lp -* because the indices of B 
are all odd integers and B ^ Lp (3 V). If the superlat- 
tice Lp is assumed, the re-scaled indices M _1 [555111] = 

[333111] are again all odd integers, so thatB £ L i p 1) {3 V) 
is again guaranteed. Yet further re-indexing of B leads 
to M" 2 [555111] = [222000], i.e., all even integers, and 
thence B and V would become equivalent through trans- 

~ (2) 

lation if the superlattice ordering were described by L p . 
Therefore, if the superlattice ordering in the approximant 
had originated from the same physical mechanism as that 
in the quasicrystal, it would have had body-centered cu- 
bic translation symmetry. Perhaps, the superlattice or- 
dering of the approximant can be understood as a kind of 
lock-in phenomenon facilitated by the discreteness of the 
atomic surfaces, while for the quasicrystal more subtle 
competitions between the boundaries of atomic surfaces 
need to be considered 2 ^. 



B. Anti-phase boundaries 

The atomic decorations of the eight half cubes that 
comprise the unit cell of the present approximant are 
mutually congruent in the sense that one of them can be 
obtained from another through a 2-fold rotation, a mirror 
reflection or the inversion. This is a strict consequence of 
the space group Pet3; in particular, two of the half cubes 
that share a square face are always connected through a 
glide operation of the space group Pa3. Although it is 
obvious that a half cube does not strictly fall on top of 
another through a half-way translation, one observes a 
partial overlap in the atomic positions between the two. 
It turns out that, from the electron diffraction pattern 
(Figure [3]) as well as the single crystal X-ray diffraction 
data, the superlattice reflections (whose Miller indices in- 
clude at least one odd number) are significantly weaker 
than the main reflections (whose Miller indices are all 
even). This indeed justifies the use of the term "superlat- 
tice ordering" in describing the present structure though 
topological rearrangements are involved in doubling the 
translational symmetry, as opposed to a more common 



type of superlattice orderings associated with density or 
displacive modulations. 

In general, the structure may not be perfectly ideal 
in realistic materials. For the present alloy sample, in 
particular, the superlattice ordering can possibly be de- 
graded to some extent due to inequilibrium processes that 
may have taken place while the material was synthesized. 
As inspired by the detailed description of the crystal 
structure (Section ^IV|) . we expect that the configuration 
of clusters within each half cube is rather strongly con- 
strained by the local packing rules, whereas the configu- 
ration of half cubes is constrained relatively loosely. In 
this respect, we find that an anti-phase boundary could 
easily be formed by introducing only a small number of 
defects in the arrangement of clusters along one of the 
glide planes^ 2 -. Hence, the material could have exhibited 
a non-negligible degradation of the superlattice ordering 
if it were obtained through a non-optimal cooling process. 
Fortunately, our sample preparation procedure (Section 
CTlf yielded a high-quality approximant crystal, in which 
only minute structural defects are included; hence, no 
special care was necessary while carrying out the struc- 
ture refinement. 

If the sample preparation were cruder, the superlattice 
ordering could be violated to a significant extent. Then, 
the satellite Bragg reflections associated with the super- 
lattice ordering could be seriously damped, possibly mak- 
ing them undetectable. Such a circumstance may have 
taken place when Sugiyama and coworkers carried out 
a structure analysis of a closely related AlyoPd23Mn6Si 
compound 6 . These authors reported that their crystal 
was a cubic 2/1 approximant with a lattice constant 
of 20.2lA and that the space group was Pmi. Note 
that the lattice constant is half times as small as that of 
the present approximant. Although the reported crystal 
structure involving a giant icosahedral cluster (diameter, 
> 20A) has never indicated a reasonable connection to 
any other crystals in closely related alloys, the knowl- 
edge was transported to the modeling of the icosahedral 
quasicrystal*^. From our revised view, the halved lattice 
constant is likely to have resulted from failing to detect 
the superlattice reflections that suffered from a serious 
damping due to anti-phase boundaries. Hence, the pre- 
vious structural solution 6 should be re-examined in view 
of that it could correspond to the average charge den- 
sity map of the eight translates {p(x + s)} of the true 
charge density map p(x), where s runs over the eight in- 
dependent translation vectors 0, R' l5 R 2 , R 3 , R^ + R 2 , 
R 2 + R' 3 , R' 3 + Ri andRi + R 2 + R 3 . 

Let us take the superposition of the corresponding 
eight translates of the true crystal structure (Subsection 
§ §IVB|) . Here, several atomic sites can fall onto the same 
position with a maximal degeneracy of 8. We find that 
cluster like arrangements are formed at the origin as well 
as at the body center of the half-sized unit cell for the su- 
perposed structure. Figure [T2l illustrates the shell struc- 
ture of the 'fake' cluster 3 at the body center, where 
most of the shells exhibit the icosahedral symmetry in 
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FIG. 12. The 'fake' cluster 3 that is observed at the body 
center of the halved unit cell in the superposition of the eight 
translated versions of the true structure. Each of the shells 
has at least the point symmetry m3 (TJ,) when the contents 
of the sites are taken into consideration. Members of a shell 
with different contents are distinguished by colors: lg (light 
gray), dg (dark gray) and bk (black). 




the geometrical arrangement except the shells 3-9 and 3- 
12. This agrees remarkably well with the shell structures 
of the primary (i.e., giant) cluster reported by Sugiyama 
et aL— . Another 'fake' cluster 03 with icosahedral ar- 
rangement is also found at the vertex of the unit cell; see 
Figure [T3] Importantly, the site degeneracies of a single 
shell in 3 or QJ generally break the icosahedral symmetry 
(Figures 12, 13 and Table IVT|) . where the minimal point 
group is m3 (Th). The violation of the icosahedral sym- 
metry in 2J may also account for the appearance of the 
secondary cluster reported by Sugiyama et o/.— . There- 
fore, the crystal structure reported by Sugiyama et al& is 
most likely to be an artifact in which the true structure 
is folded into half-sized unit cell. 



C. Atomic packing 

As described in Subsection ^IVBl the closest inter- 
atomic distances in the idealized construction are given 
by b' and d and are parallel to two- and three-fold axes, 
respectively. These distances are 1/r 2 times smaller than 
b- and c-linkages for the cluster packing, so that the ratio 
between the two satisfies c'/b' — c/b (= y3/2). Then, it 
is natural that atoms would be locally arranged to sit on 
the vertices of the four basic shapes, A', B', C and D', 



FIG. 13. The 'fake' cluster 2J that is observed at the origin 
in the superposition of the eight translated versions of the 
true structure. Each of the shells has at least the point sym- 
metry m3 (Th) when the contents of the sites are taken into 
consideration. Members of a shell with different contents are 
distinguished by colors: lg (light gray), dg (dark gray) and bk 
(black). 
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which are the miniature versions (scale, xl/r 2 ) of A-, 
B-, C- and D-cells, respectively. Yet, detailed inspection 
reveals several essential differences between the arrange- 
ment of atoms and that of clusters; that is, the former is 
not precisely a miniaturized version of the latter. 

An important feature of the atomic packing is that 
there exists an icosahedral configuration located at the 
center of each B-clustcr, consisting of the central site 
(type, (Bo)) and the 12 surrounding sites (type, (B 5 ,...)) 
forming the inner icosahedral shell. Hence, it is required 
to introduce an icosahedron with a central dot as a ba- 
sic shape if a tiling description of the atomic packing is 
sought. We denote this shape I' (Figure [TJ] left). 

Remember also that there exist two splitting sites 
(type, (M 2 )) associated with each type I b-linkage that 
connects two M-clusters. Although the two sites cannot 
be occupied simultaneously, an Al atom is most likely 
jumping back and forth between these sites at a typical 
annealing temperature of 850° C, where both the splitting 
sites play an inevitable role in maintaining and stabilizing 
the structure. The two splitting sites are always located 
at the two tips of a flat hexagonal bipyramid, which cor- 
responds to the overlap of two adjacent M-clusters. We 
denote this shape H' (Figure [HI right). The two tips 
of H' are the very splitting sites of the type (M2), while 
the six vertices of the base hexagon include two atomic 
sites that are symbolized as (M 5 ,M 5 ,...) and four that are 
symbolized as (M2,M2,...). 

The idealized atomic arrangement can now be repre- 
sented as an arrangement of the six basic geometrical 
shapes A', B', C, D', I' and IT. Still, it may not be 
called a tiling because a part of the space is left blank 
in this description. The latter circumstance arises when 
the cluster packing involves any number of D-cells. Re- 
member that a D-cell has the largest volume among the 
four canonical cells and clusters are free from interpene- 
tration at its central region. This causes that the inner 
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TABLE VI. The content of each site in the fake clusters 3 and 
QJ, which are depicted in Figures 12 and 13, respectively. The 
shells are numbered in increasing order of radius, where the 
central site of QJ is numbered 0. 

shell # radius content 

[A] 



(cluster 3) 



8 
9 

10 



11 



12 
13 



14 



15 



16 



1.738 
2.813 
4.144 
4.551 
5.550 
6.330 



2(B 5 ) 

2(B 3 M 5 ), 4<B§M 5 > 
2(M|), 6(B 5 M|) 
2(B ) 
2<B 3 M 3 ) 

lg: (B 5 M|), (B 5 M 2 ), 4(M 2 ) 
bk: 2(B 5 M 2 ), 4(M 2 ) 



6.705 lg: 6<BiM^>, 2<B 3 M|> 
bk: 3<B|M|), 5(B 3 Ml) 
2(B 5 M 2 ), 6(M ) 
<M 2 ) 



7.364 
7.935 
8.237 



dg: (B^M 5 >, <B^>, 6<B 3 M 3 ) 



8.782 



9.526 
10.001 



lg: 
bk 
(lg: 
lg: 

bk 



<B§>, 5(B 3 M 3 ) 
2(B|M 5 ), 6(B 3 M 3 ) 
4(B 5 M|), 4(B 5 M 2 ) 
6(B 5 Mi), 2(B 5 M 2 ) 
3<B 5 Ml), (Mi), 2(B 5 M 2 ), 2(M 2 ) 



2(B 3 M 3 ) 



10.242 



10.685 



10.849 



dg: 
lg: 
bk: 
dg: 
lg: 
bk: 
dg: 
lg: 
bk: 
lg: 
bk: 



2(Ml), 4(M 2 ) 
(M 2 ) 

(Ml), 3(M 2 ) 
4<B§M 5 >, 4{B|Ml> 
(B 3 M 3 ), (B§), 4{B|M 5 ), 2(B|M|) 
<B 3 M 3 ), 3<BiM 5 ), 4<B§M|> 
2(B 5 M 2 ), 6<B 5 Ml) 



3 5 Ml) 



(M ), 
8<B ) 
5(B ) 



:B 5 M 2 ), 6(B 5 M^) 



(cluster QJ) 



0.000 8(M ) 
2.561 lg: 2(B 3 M 3 ) 
bk: 5<B 3 M 3 ) 
2.957 lg: 4(M 2 ) 
bk: <M 2 ) 
8<BlMi) 



4.551 
4.785 

6.604 



.705 



lg: 4{B 5 Ml) : 



bk 
dg: 
lg: 
bk 
lg: 
bk 



7<B 5 M|) 
(B3M5}, 
4<BiM 5 ), 
(B|M 5 ), 
2(B ) 
5(B ) 



4(M 2 ) 
, (M 2 ) 
5(B 3 M 3 ) 
2<B 3 M 3 ) 
2<B 3 M 3 ) 



shell atoms of an M-cluster located at one of the vertices 
are sparsely distributed toward the center of the D-cell, 
thereby providing a large unfilled space or a void inside 
the outer shell of the M-cluster. Two relevant kinds of 
void are illustrated in Figure [TS] for the type I and type 
II D-cells. These voids can be interconnected to form 



FIG. 14. The two basic tiles, I' (left) and H' (right), that 
are necessary to provide a tiling description for the atomic 
positions. The double bars represent the b'-linkages, while 
the thick bars the c'-linkages between atoms. There is an 
atomic site for Pd at the center of the tile I'. The surfaces 
of the tiles are composed of X'- and Y'-faces, which are the 
miniature versions of the X- and Y-faces, respectively. + 
or — signs shown on the front faces indicate which of the 
symmetrically distinguishable sides are facing outward. Note 
that the point symmetry of I' and H' are 532/m (7h) and 
mmm (-D 2 /i), respectively. 




a larger void or a channel if D-cells adjoin together. In 
the refinement of the present approximant, it is observed 
that a surrounding atom of such a void can have a plit- 
ting position within the void (Pdl57' for Figure[15ja) and 
A186' for Figure Hl^ b)) or exhibit a significant displace- 
ment toward the center of the void (All 74' for Figure 
IT5T b)). These facts are physically consistent with the 
present geometrical description. 



VI. CONCLUDING REMARKS 

A simple heuristic search in the quaternary Al-Pd-Cr- 
Fe system has led to the first discovery of a highly stable 
cubic approximant phase to Al-based F-type icosahedral 
quasicrystals i-Al7oPd 2 oTMio (TM=transition metals). 
Single crystals of the approximant a few hundred /urn's 
in diameter have been grown with slow cooling, suggest- 
ing that the approximant phase is formed congruently 
from the melt. Using single crystal X-ray diffraction, an 
ab initio structure analysis has been carried out success- 
fully, revealing a crystal symmetry of Pa3 with a lattice 
constant of 40.5 A. All the atoms in the crystal belong to 
two kinds of clusters which in the literature are referred 
to as the Mackay type and the Bergman type clusters, 
whereas there is no need for glue atoms. 

The centers of the clusters are located at the vertices of 
a CCT, designated as a 2x2x2 superstructure of the 3/2- 
packing^, in which the parity of each vertex determines 
unambiguously the kind of the relevant cluster. Impor- 
tantly, the present structure offers a set of universal rules 
for the local packing of clusters. The rules are applica- 
ble to arbitrary CCT's, whereby a number of hypothet- 
ical approximants can be constructed. Some of these 
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FIG. 15. Two kinds of voids associated with D-cells. (a) A 
void within the outer shell of an M-cluster located at a corner 
D a of a type I D-cell. The edges of the D-cell are represented 
as black bars. The spheres are the atomic positions predicted 
by the present geometrical rules, where labels for some atomic 
sites are presented in the figure, (b) Similar to (a) but the void 
is associated with a corner Dt of a type II D-cell. A splitting 
position for an atomic site is indicated with an arrow and a 
dotted circle (labeled 's'). A possible atomic displacement is 
indicated with an arrow (labeled 'd'). In the refinement, the 
two kinds of splitting positions as in (a) and (b) are included 
as Pdl57' and A186', respectively. 




hypothetical structures could be synthesized in related 
alloys. Indeed, the crystal structure that was previously 
reported as a cubic 1/1 approximant by Sugiyama et al£ 
could be better understood as an artifact of folding a 
hypothetical structure constructed from a 2x2x2 super- 
structure of the 2/1 packing^ into a half-sized unit cell. 
Here a similar argument to that presented in Subsection 
§ WBI is assumed. A search for stable approximants with 
various degrees of rational approximation to r within re- 
lated alloy systems is an important future challenge. 

In actual fact, the proposed two cluster units, M and B, 
have been already used (with a somewhat different way 
of handling the inner shell of M-cluster) to compose a 
hypothetical model of F-type icosahedral quasicrystals, 
e.g., ^-Al-Pd-M n 23 ' 24 . In this case, the clusters are as- 



sumed to be centered at the vertices of a random tiling 
composed of two kinds of rhombohcdra called Ammann's 
rhombohcdra. Accordingly individual clusters are con- 
nected with their adjacent neighbors through the edges 
(in five-fold directions) or the short diagonals (in two-fold 
directions) of the rhombic faces, whereas clusters at the 
two ends of the short body diagonal (along a three-fold 
direction) of each obtuse rhombohedron are forced to be 
interpenetrated heavily with each other. Our findings 
proved to be against this model, since no pair of clusters 
connected through a five-fold linkage is observed in the 
reconstructed crystal structure. 

The superlattice ordering in the icosahedral quasicrys- 
tal, i-Al-Pd-Mn^ has once attracted much attention, and 
rich features associated with this phenomenon have been 
published 2 — ~—. However, no clear account on its phys- 
ical origin has been provided, although it seems to be 
closely related to the binary cluster composition for the 
F-type icosahedral systems. In the expectation that the 
superlattice ordering in the present approximant may be 
related to that in the quasicrystal, possible P-type su- 
perlattice orderings in the six-dimensional F-type hyper- 
cubic lattice have been analyzed in Subsection A[ 
where it has turned out that the two phenomena cannot 
be associated with the same six-dimensional superlattice. 
Perhaps, a kind of lock-in transition facilitated by the dis- 
creteness of the atomic surfaces is more likely to be the 
origin of the superlattice ordering in the approximant, 
whereas in the quasicrystal more subtle competitions be- 
tween the boundaries of atomic surfaces could be respon- 
sible for it. 

One of the most important outcomes of this work is 
the finding that the atomic arrangement is perfectly de- 
scribed using a set of geometrical rules for cluster packing 
defined on a CCT. Local arrangements of atomic posi- 
tions thus determined can be described generally as the 
vertices of six basic polyhedra A', B', C, D', I' and H', 
which are arranged under the face-to-face matching con- 
straint. Importantly, the latter construction leaves the 
presence of a void inside the outer shell of each M-cluster 
that is located at a vertex of a D-cell. The present struc- 
ture refinement has indicated that some of the surround- 
ing atoms of the voids have a secondary (i.e., splitting) 
position within the void or exhibit a significant displace- 
ment toward the center of the void. 

Altogether, the present results have greatly improved 
our understandings on the local atomic structures of F- 
type icosahedral quasicrystals and their approximants. 
Clearly, the present local rules for cluster packing main- 
tain a certain level of universality, and a variety of com- 
plex structures that are closely related to Al-based F- 
type icosahedral quasicrystals are expected to be de- 
scribed within the same framework. However, when the 
structure of the quasicrystal is at stake, there remains a 
crucial issue of whether the canonical cells can be used 
to construct a quasiperiodic tiling with the icosahedral 
symmetry. The peculiar shapes of the canonical cells 
are known to impose non-local geometrical constraints, 
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making it difficult to determine their non-periodic align- 
ments throughout the whole space. It is an important 
challenge in the theoretical crystallography of quasicrys- 
tals to prove if an icosahedral quasiperiodic tiling can be 
constructed using the canonical cells. 



Appendix A: Derivation of the statistics formulae 
for the 16 types of atomic sites 



Note that each class of atomic sites is associated with 
one or more type(s) of geometric objects (see Figure [5] 
and Table IIVI) , and a careful inspection of the decorated 
objects proves the following relationships: 



n(nodei) 


= f 






(Al) 




= Oh 


-m)/20 




(A2) 




= (n- 


h 3r + s + It •+ 


■3u)/12 


(A3) 




= (<H 


-2k + l + 2p + 


■3g)/30, 


(A4) 


n(noden) 


= i 






(A5) 




= CH 


-l+p)/12 




(A6) 




= (m 


+ n + 2o + 2s 


+ 2i + u)/20, 


(A7) 


n(bi) 


= 9/2 






(A8) 




= (3r 


+ < + 3u)/2 




(A9) 




= (AH 


-p + 3 ? )/4, 




(A10) 


n(bn) 


= ( 4 


- « + 0/2, 




(All) 


ra(c) 


= m 






(A12) 




= (l + 


2p)/3 




(A13) 




= (n- 


h 2s + At + 3zt)/3, 


(A14) 


rc(Xj) 


= P 






(A15) 




= 2t-f 


- 3u, 




(A16) 


n(X n ) 


= s + 


2t, 




(A17) 


n(Y x ) 


= q 






(A18) 




= r + 


u, 




(A19) 


n(Ai) 


= t, 






(A20) 


n(d) 


= u, 






(A21) 



where the variables /, g, u represent the number fre- 
quencies for the 16 types of atomic sites; see Table IIVI 
for the definitions. Some of these equations turn out to 
be redundant in the sense they are not necessarily in- 
dependent with each other. By solving these equations, 
Eas.(|32jl - (471) of Subsection § JVC] can be readily ob- 
tained. 



Appendix B: Refined parameters 

Among the 4728 atomic positions contained in a sin- 
gle unit cell, 4680 are reproduced (within small devia- 
tions) by replicating the cluster templates at the vertices 
of the CCT under the geometrical rules fully described in 
Subsection § §1VBI The remaining 48 positions are irreg- 
ular splitting sites (A186' and Pdl57') introduced inside 
voids associated with M-clusters centered at Fe/Pd82 and 
Cr/A198. The inclusion of these splitting sites has re- 
duced R W (F) by more than 1%. Bear in mind that no 
glue atom is included in the structure. 

The existence of voids causes some additional irregu- 
larities for the refinement: (1) The splitting site Pdl57' 
appears to contradict with the site Al/Crl04 which lies 
nearby. Therefore, the total s.o.f. of Al/Crl04 is set 
to be equal to that of Pdl57, so that it complements 
that of Pdl57'. (2) In addition, the atomic site All 74' 
shifted toward the center of a void lies too close to A1183, 
whereas the distance between All 74 and All 74', which 
are originally splitting sites for a single Al atom, is large 
enough to accommodate atoms at both of them. There- 
fore, s.o.f.'s for the four positions A1174, A1174', AU83, 
A1183' are refined so that their sum is fixed to 2. These 
are exceptional measured that we have taken to achieve 
reasonable convergence of the refinement. 

The number of independent atomic positions cannot 
be reduced below 204. As an effort to suppress the num- 
ber of independent parameters, the following three con- 
straints on the isotropic atomic displacement parameters 
(IADP's) are introduced. (1) Within the inner shell of 
each B-cluster, the IADP's are set to be identical for the 
same site type. See, for example, the IADP's of A112, 
A123 in Table ED] the IADP's for the site classes 
(B 5 ,M 2 ,M 2 ), (B 5 ,M 2 ) and (B 5 ) are 0.015, 0.01 and 0.006, 
respectively. (2) The IADP's for the two splitting sites 
associated with each b- linkage between M-clusters are set 
to be identical; e.g., the IADP's for All 73 and A1173' are 
both 0.029. (3) In the final stage of the refinement, it 
was found that the IADP and the site occupation factor 
(s.o.f.) are strongly correlated for Fe/Pdl54. If both the 
parameters are refined, the IADP converges to a negative 
value while the Pd concentration reduces by more than 
10 %. Hence, the IADP of this atomic site is fixed to be 
0.011. 

In Table IVIII it is observed that IADP's for some 
atomic positions in the inner shell of M-clusters are ex- 
ceptionally large (e.g., A196, A197, A100 and Al/Crl04). 
This might be a due to possible disorders in the M3 shells. 
A full characterization of such disorders may require a de- 
tailed analysis of charge density distributions for the M3 
shells. However, the quality of our intensity data may 
not be sufficient for performing such an analysis. 



TABLE VII: Table for atomic parameters for the refined structure, s.o.f. 
gives the site occupation factor for each atomic species. U is the IADP. 
A gives the distance between the refined atomic position and the ideal 
position. 



atom 


Wyckoff 


s.o.f. 


x/aiat 


y/aiat 


z/ai at 


U 


site type 




Pdl 


8c 


1.00 


0.09485(8) 


0.09485(8) 


0.09485(8) 


0.003 


/r>U)\ 


0.001 


A12 


24d 


1.00 


0.1308(3) 


0.0944(3) 


0.0381(3) 


0.005 


/Ti(l) TV T ( 3 ) \ 


0.002 


A13 


24d 


1.00 


0.0585(4) 


0.0946(4) 


0.0386(3) 


0.005 


/i-i(l) at(1) tv t ( 3 ) \ 

(B 5 ; ,M^ ; ,M^ ; ) 


0.002 


A14 


24d 


1.00 


0.1516(4) 


0.1300(3) 


0.0949(4) 


0.005 


/Ti(l) TV X ( 5 ) 1 f(5)\ 

(B K B ; ,M^ 'M 2 } 


0.004 


A15 


24d 


1.00 


0.0948(4) 


0.1515(4) 


0.0597(4) 


0.005 


/Ti(l) TV T ( 3 ) 1 f(5)\ 


0.003 


Pd6 


8c 


1.00 


0.40400(9) 


0.40400(9) 


0.40400(9) 


0.008 


/t>(2)\ 


0.001 


A17 


24d 


1.00 


0.4612(4) 


0.4387(4) 


0.4051(4) 


0.012 


(B 5 'Ml 'M),'} 


0.003 


A18 


24d 


1.00 


0.4043(4) 


0.4619(3) 


0.3680(3) 


0.002 




0.002 


A19 


24d 


1.00 


0.3698(4) 


0.4025(4) 


r\ o A f o / A\ 

0.3458(4) 


0.012 


(B 5 ; ,M^ 'M 2 ) 


0.003 


A110 


24d 


1.00 


0.4389(4) 


0.4046(4) 


0.3487(4) 


0.012 


(B 5 ; ,M^ 'M 2 ) 


0.004 


Pdll 


24d 


1.00 


0.25036(9) 


0.3449(1) 


0.19066(9) 


0.006 


/t>(3)\ 
< B o ') 


0.001 


A112 


24d 


1.00 


0.2868(4) 


0.3453(4) 


0.1326(4) 


0.015 


/TT>(3) if(3) at(6)\ 


0.001 


A113 


24d 


1.00 


0.3459(4) 


0.2472(4) 


0.2142(5) 


0.015 


/r>(3) ,,(6) ,,(T)i 

(B^ 'Ml 'Ml') 


0.003 


A114 


24d 


1.00 


0.2503(3) 


0.4026(4) 


0.2255(4) 


0.010 


/TTi(3) TV T ( 7 ) \ 

<B^ \M^) 


0.003 


A115 


24d 


1.00 


0.1915(4) 


0.3809(4) 


0.1905(4) 


0.010 


/TTi (3) TV t(6) \ 

(B5 ; ,M^ ; ) 


0.001 


A116 


24d 


1.00 


0.3074(4) 


0.3101(4) 


0.1894(4) 


0.015 


/t>(3) ,,(6) ,,(r)\ 
(B^ ; ,M^ 'Ml ) 


0.003 


A117 


24d 


1.00 


0.2508(3) 


0.2898(4) 


0.2264(3) 


0.006 


(B5 > 


0.004 


A118 


24d 


1.00 


0.3074(4) 


0.3804(4) 


0.1899(4) 


0.015 


/1-> (3) Tv T ( 3 ) 1 r(7)\ 

(Bg ',Mj ,M'2 ) 


0.002 


A119 


24d 


1.00 


0.2160(4) 


0.3474(4) 


0.1333(4) 


0.010 


<B^ 'Ml') 


0.003 


A120 


24d 


1.00 


0.2490(4) 


0.4013(4) 


0.1564(4) 


0.010 


/t~» (3) n /r ( 3 ) \ 

(B 5 >,M { 2 >) 


0.004 


A121 


24d 


1.00 


0.2861(4) 


0.3468(4) 


0.2471(4) 


0.015 


/ t-\ (3) ti r(tt ti r(tt\ 

(B y 5 >,M y 2 >,M 2 ') 


0.003 


A122 


24d 


1.00 


0.1935(4) 


0.3102(4) 


0.1902(4) 


0.015 


<B£ ; ,M^ ; ,M^ ; ) 


0.003 


A123 


24d 


1.00 


0.2492(4) 


0.2876(5) 


0.1557(4) 


0.015 


(B^ ; ,M^ ; ,M^ ; ) 


0.002 


Pd24 


24d 


1.00 


0.34558(9) 


0.49947(11) 


0.2496(1) 


0.008 




0.001 


A125 


24d 


1.00 


0.2514(4) 


0.2123(4) 


0.0366(4) 


0.010 


/Tr>(4) TV j(5) T. t(6) \ 

(B^ ; ,M^ ; ,M^ ; ) 


0.002 


A126 


24d 


1.00 


0.4019(5) 


0.4640(4) 


0.2491(4) 


0.010 


/TTi(4) TV T (4) tv t(7) \ 


0.003 


A127 


24d 


1.00 


0.3793(4) 


0.4985(4) 


0.3065(4) 


0.012 


/TTi(4) TV t(4)\ 

<B^ j M ) 


0.004 


A128 


24d 


1.00 


0.2881(4) 


0.4654(4) 


0.2491(4) 


0.012 


/t->(4) a/t(7)\ 
< B 5 ; ) 


0.003 


A129 


24d 


1.00 


0.1930(4) 


0.1181(4) 


-0.0004(4) 


0.012 


/Ti(4) tvt(3)\ 

(B^ ; ,M^ ; ) 


0.002 


A130 


24d 


1.00 


0.3461(4) 


0.4437(4) 


0.2840(4) 


0.010 


/n(4) tv t(7) a t(7) \ 

(B^ ; ,M^ ; ,M^ j ) 


0.004 


A131 


24d 


1.00 


0.3457(4) 


0.4426(4) 


0.2158(4) 


0.010 


/n(4) ,,(3) ,,(7), 

(B^ ; ,M^ ; ,M^ ; ) 


0.003 


A132 


24d 


1.00 


0.2834(4) 


0.1540(4) 


0.0569(4) 


0.010 


/Ti(4) TV »-(5) TV t(6) \ 

(B^ >,M 2 'Ml') 


0.004 


A133 


24d 


1.00 


0.2132(4) 


0.1562(4) 


0.0572(4) 


0.010 


/i-i(4) tv r (5) tv r (5) \ 


0.003 


A134 


24d 


1.00 


0.1921(4) 


0.1897(4) 


0.0007(4) 


0.010 


(B^ >,M 2 'Ml') 


0.002 


A135 


24d 


1.00 


0.3056(4) 


0.1904(4) 


-0.0001(4) 


0.010 


/Ti(4) tv t(6) tv tr(7) \ 

{B\ 'Ml 'M 2 ') 


0.003 


A136 


24d 


1.00 


0.2489(4) 


0.0980(5) 


r\ r\ f rv / a \ 

0.0352(4) 


0.010 


/td(4) ,j(4) „(5), 


0.003 


Pd37 


24d 


1.00 


0.40403(9) 


0.09507(9) 


0.09482(9) 


0.006 


(B^ ; ) 


0.001 


A138 


24d 


1.00 


0.4391(4) 


0.0940(4) 


0.0384(4) 


0.007 


(Bg ' Ml Ml ) 


0.003 


A139 


24d 


1.00 


0.1299(4) 


0.4040(4) 


0.0383(4) 


0.008 


/t>(8) ,,(4)i 
(Bg ; ,M^ ; ) 


0.003 


A140 


24d 


1.00 


0.3692(4) 


0.0969(4) 


0.0380(3) 


0.008 


(B^M { 2 } ) 


0.002 


A141 


24d 


1.00 


0.0595(4) 


0.4051(4) 


0.0375(4) 


0.007 


(B 5 5 \M^,M( 4) > 


0.001 


A142 


24d 


1.00 


0.0944(4) 


0.4614(4) 


0.0600(4) 


0.007 


(B 5 5 ),M( 2) ,M( 4) > 


0.003 


A143 


24d 


1.00 


0.4616(4) 


0.1302(4) 


0.0957(4) 


0.008 


(BfM 4 ') 


0.003 


A144 


24d 


1.00 


0.4043(4) 


0.1518(4) 


0.0607(4) 


0.007 


(B 5 5 ),M( 4) ,M( 6) ) 


0.003 


A145 


24d 


1.00 


0.1523(4) 


0.3699(4) 


0.0951(4) 


0.008 


<Bf\M< 5) > 


0.003 


A146 


24d 


1.00 


0.3477(4) 


0.1309(4) 


0.0954(4) 


0.007 


(b 5 5) ,m( 5) ,m( 6) ) 


0.002 


A147 


24d 


1.00 


0.4047(4) 


0.1521(4) 


0.1314(4) 


0.008 


(Bf ,Mf) 


0.002 


A148 


24d 


1.00 


0.1520(4) 


0.4391(4) 


0.0946(4) 


0.008 


(B ( 6 5) M ( 2 4) ) 


0.003 
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A149 


24d 


1.00 


0.0951(4) 


0.3472(4) 


0.0605(4) 


0.007 


(Bf ,H4 4 \MW ) 


0.002 


Pd50 


24d 


1.00 


0.5000(1) 


0.25064(9) 


0.15386(9) 


0.004 


<B< 6) ) 


0.001 


A151 


24d 


1.00 


0.3062(4) 


0.3107(4) 


0.0003(4) 


0.009 


(B^M^M^) 


0.003 


A152 


24d 


1.00 


0.4989(4) 


0.1934(4) 


0.1882(4) 


0.009 


<B< 6) ,M< 5) ) 


0.004 


A153 


24d 


1.00 


0.4638(4) 


0.2495(4) 


0.0970(4) 


0.009 


(Bf,M( 4) M 6) > 


0.002 


A154 


24d 


1.00 


0.4985(4) 


0.3065(4) 


0.1190(4) 


0.009 


<Bf M 3) > 


0.003 


A155 


24d 


1.00 


0.5359(4) 


0.2490(4) 


0.2103(4) 


0.009 


<B< 6) ,M< s) ,m( 6) > 


0.003 


A156 


24d 


1.00 


0.4998(4) 


0.1930(4) 


0.1203(4) 


0.009 


<B< 6) ,M< 4) ) 


0.003 


A157 


24d 


1.00 


0.4436(4) 


0.2836(4) 


0.1542(4) 


0.009 


(Bf ,MfM 7) > 


0.004 


A158 


24d 


1.00 


0.4417(4) 


0.2141(4) 


0.1543(4) 


0.009 


(Bf ,Mf) 


0.001 


A159 


24d 


1.00 


0.2151(4) 


0.3465(4) 


0.0582(4) 


0.009 


(Bf ,M ( 2 } ) 


0.002 


A160 


24d 


1.00 


0.2841(4) 


0.3448(4) 


0.0568(4) 


0.009 




0.003 


A161 


24d 


1.00 


0.2502(4) 


0.4030(4) 


0.0342(4) 


0.009 


<B< 6) ,M< 3) ,M< 4) ) 


0.003 


A162 


24d 


1.00 


0.4650(4) 


0.2496(4) 


0.2105(4) 


0.009 


<Bf ,M< 7) > 


0.003 


Pd63 


24d 


1.00 


0.49994(9) 


0.44058(9) 


0.15447(8) 


0.007 


<B^ 7) ) 


0.000 


A164 


24d 


1.00 


0.3822(4) 


0.3091(4) 


0.0010(4) 


0.012 


<Bf\M( 6 \lv4 7) > 


0.001 


A165 


24d 


1.00 


0.4633(4) 


0.4390(4) 


0.2120(4) 


0.012 


<b( 7) ,m< 4) ,m< 7) > 


0.003 


A166 


24d 


1.00 


0.5340(4) 


0.4421(4) 


0.2129(4) 


0.012 




0.003 


A167 


24d 


1.00 


0.4425(4) 


0.4024(4) 


0.0353(4) 


0.012 


(b«,m«,m( 3) ) 


0.003 


A168 


24d 


1.00 


0.4756(3) 


0.3449(3) 


0.0579(4) 


0.008 


<b«m( 3) > 


0.002 


A169 


24d 


1.00 


0.4431(4) 


0.4768(3) 


0.1556(3) 


0.008 


<B< 7) ,M< 3) > 


0.002 


A170 


24d 


1.00 


0.4986(4) 


0.3836(3) 


0.1198(3) 


0.008 


<Bf\]v4 3) ) 


0.003 


A171 


24d 


1.00 


0.5003(4) 


0.4976(4) 


0.1184(3) 


0.012 




0.002 


A172 


24d 


1.00 


0.4974(4) 


0.3108(3) 


0.0004(5) 


0.012 


<Bf\]v4 3) ,]v4 4) } 


0.003 


A173 


24d 


1.00 


0.4439(4) 


0.4055(4) 


0.1546(4) 


0.012 


<b«m( 3 \m«> 


0.003 


A174 


24d 


1.00 


0.4051(5) 


0.3453(4) 


0.0567(4) 


0.012 


(b«,m( 3) ,m< 6) ) 


0.002 


A175 


24d 


1.00 


0.4653(4) 


0.4400(4) 


0.0964(4) 


0.012 


<B«M«Mf ) 
<M«> 


0.002 


Cr/A176 


4a 


0.76/0.24 











0.015 


0.000 


Cr/A177 


24d 


0.88/0.12 


0.47737(18) 


0.44032(18) 


-0.0002(2) 


0.007 


<B«M«> 


0.001 


Cr/A178 


8c 


0.96/0.04 


0.03635(18) 


0.03635(18) 


0.03635(18) 


0.006 


(B«,M«) 


0.000 


Cr/A179 


4b 


0.54/0.46 


0.50000 








0.007 


<Mf> 


0.000 


Fc/Pd80 


8c 


0.75/0.25 


0.46330(13) 


0.46330(13) 


0.46330(13) 


0.002 


<B| 2) ,Mf> 


0.000 


Fe/Pd81 


24d 


0.76/0.24 


0.46332(15) 


0.03661(15) 


0.03672(15) 


0.013 


<B| 5) ,Mf> 


0.000 


Fe/Pd82 


24d 


0.89/0.11 


0.34520(19) 


0.44053(19) 


0.0954(2) 


0.017 


<M< 3) > 


0.001 


A183 


24d 


1.00 


0.3114(5) 


0.4078(6) 


0.0598(5) 


0.048 


<B< 6 \M< 3 >> 


0.004 


A184 


24d 


1.00 


0.3817(6) 


0.4086(7) 


0.0601(7) 


0.081 


<B«Mf> 


0.004 


A185 


24d 


1.00 


0.3458(5) 


0.4624(5) 


0.0387(5) 


0.045 


<B«Mf> 


0.003 


A186 


24d 


0.45 


0.3078(8) 


0.4079(8) 


0.1284(8) 


0.008 


(Bf,Mf) 


0.005 


A186' 


24d 


0.55 


0.2911(6) 


0.4358(6) 


0.1186(6) 


0.008 


(Bf,M< 3) ) 


0.038 


A187 


24d 


1.00 


0.3447(5) 


0.4619(5) 


0.1520(5) 


0.025 


(B| 4) ,Mf) 


0.003 


A188 


24d 


1.00 


0.4008(6) 


0.4398(6) 


0.1181(6) 


0.077 


(B«,Mf) 


0.004 


A189 


24d 


1.00 


0.3692(6) 


0.4957(6) 


0.0937(6) 


0.075 


<B«Mf> 


0.005 


Fe/Pd90 


24d 


0.91/0.09 


0.4990(2) 


0.19166(13) 


0.0003(2) 


0.012 


<M< 4 >> 


0.001 


A191 


24d 


1.00 


0.2894(5) 


0.0563(5) 


0.0003(5) 


0.041 


(B ( 3 4) ,M ( 3 4) ) 


0.004 


A192 


24d 


1.00 


0.2136(5) 


0.4439(6) 


0.0026(6) 


0.054 


{B ( 3 &) ,M ( 3 4) ) 


0.004 


A193 


24d 


1.00 


0.3446(4) 


0.0345(4) 


0.0338(4) 


0.007 


(b| 5) ,m| 4) ) 


0.003 


A194 


24d 


1.00 


0.1584(6) 


0.4661(6) 


0.0348(6) 


0.065 


<b| 5 \m< 4 >> 


0.005 


A195 


24d 


1.00 


0.4991(6) 


0.4817(5) 


0.2521(5) 


0.048 


(Bi 7) ,M< 4) ) 


0.005 


A196 


24d 


1.00 


0.4670(8) 


0.4680(8) 


0.3470(8) 


0.105 


(Bf,Mi 4) ) 


0.006 


A197 


24d 


1.00 


0.469(1) 


0.153(1) 


0.032(1) 


0.152 


(B| 5) ,Mi 4 )) 


0.007 


Cr/A198 


24d 


0.80/0.20 


0.1531(2) 


0.2497(2) 


0.0947(2) 


0.001 


<M< 5) > 


0.002 
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atom 


Wyckoff 


s.o.f. 


x/ aiat 


y/a-iat 


z/ai at 


Cr/A199 


24d 


0.72/0.28 


0.1918(2) 


0.2162(2) 


0.0607(2) 


A1100 


24d 


1.00 


0.1748(8) 


0.2803(8) 


0.0502(8) 


Cr/AllOl 


24d 


0.85/0.15 


0.13310(19) 


0.1906(2) 


0.0953(2) 


Cr/A1102 


24d 


0.79/0.21 


0.13219(19) 


0.3076(2) 


0.0952(2) 


Cr/A1103 


24d 


0.77/0.23 


0.2480(2) 


0.1194(2) 


0.0962(2) 


Al/Crl04 


24d 


0.38/0.27 


0.191(1) 


0.290(1) 


0.116(1) 


Fe/Pdl05 


24d 


0.98/0.02 


0.1555(2) 


0.75056(19) 


0.0956(3) 


A1106 


24d 


1.00 


0.3088(4) 


0.2842(4) 


0.1317(4) 


Al/Crl07 


24d 


0.77/0.23 


0.3087(4) 


0.2844(4) 


0.0586(4) 


Al/Crl08 


24d 


0.83/0.17 


0.3837(4) 


0.2834(4) 


0.0608(4) 


Al/Crl09 


24d 


0.89/0.11 


0.4034(4) 


0.2477(4) 


0.1180(4) 


A1110 


24d 


1.00 


0.3087(5) 


0.2128(5) 


0.0597(5) 


Aim 


24d 


1.00 


0.3446(5) 


0.2258(5) 


0.1524(4) 


Aim 


24d 


1.00 


0.3672(6) 


0.1931(6) 


0.0976(6) 


Fell3 


24d 


1.00 


0.4047(3) 


0.3455(2) 


0.2498(2) 


A1114 


24d 


1.00 


0.4421(4) 


0.3802(4) 


0.2141(4) 


A1115 


24d 


1.00 


0.4054(4) 


0.3663(4) 


0.3083(4) 


A1116 


24d 


1.00 


0.3846(4) 


0.4039(4) 


0.2495(4) 


A1117 


24d 


1.00 


0.3464(4) 


0.3445(4) 


0.2263(4) 


A1118 


24d 


1.00 


0.4392(4) 


0.3104(4) 


0.2143(4) 


A1119 


24d 


1.00 


0.4388(5) 


0.3102(5) 


0.2841(5) 


A1120 


24d 


1.00 


0.3823(4) 


0.2882(4) 


0.2501(4) 


Pdl21 


24d 


1.00 


0.0379(1) 


0.1538(1) 


0.0334(1) 


Pdl22 


24d 


1.00 


0.15498(9) 


0.34383(9) 


0.0383(1) 


Pdl23 


24d 


1.00 


0.4998(1) 


0.40696(9) 


0.36725(9) 


Pdl24 


24d 


1.00 


0.46348(11) 


0.15417(11) 


0.1535(1) 


Pdl25 


24d 


1.00 


0.4431(1) 


0.49909(11) 


0.2163(1) 


Pdl26 


24d 


1.00 


0.43854(9) 


0.38408(9) 


0.0950(1) 


Pdl27 


24d 


1.00 


0.4613(1) 


0.3446(1) 


0.3433(1) 


Pdl28 


24d 


1.00 


0.25334(9) 


0.38382(9) 


0.0950(1) 


Pdl29 


24d 


1.00 


0.34485(11) 


0.15618(11) 


0.03799(11) 


Pdl30 


24d 


1.00 


0.21331(9) 


0.44028(9) 


0.19174(9) 


Pdl31 


24d 


1.00 


0.09127(9) 


0.36632(9) 


0.00029(11) 


Pdl32 


24d 


1.00 


0.28893(9) 


0.4407(1) 


0.18784(9) 


Pdl33 


24d 


1.00 


0.1904(1) 


0.40743(11) 


0.0554(1) 


Pdl34 


24d 


1.00 


0.19072(9) 


0.40400(9) 


0.13116(9) 


Pdl35 


24d 


1.00 


0.31200(11) 


0.28740(11) 


0.2495(1) 


Pdl36 


24d 


1.00 


0.3089(1) 


0.09278(11) 


0.0566(1) 


Pdl37 


24d 


1.00 


0.1526(1) 


0.1546(1) 


0.03373(11) 


Pdl38 


24d 


1.00 


0.2123(1) 


0.25038(9) 


0.00257(11) 


Pdl39 


24d 


1.00 


0.4390(1) 


0.3074(1) 


0.09494(11) 


Pdl40 


24d 


1.00 


0.1913(1) 


0.0961(1) 


0.0599(1) 


Pdl41 


24d 


1.00 


0.46242(11) 


0.3452(1) 


0.1567(1) 


Pd/Fel42 


24d 


0.74/0.26 


0.28543(11) 


0.24948(11) 


0.19141(12) 


Fe/Pdl43 


24d 


0.51/0.49 


0.50141(13) 


0.28711(12) 


0.05893(12) 


Pd/Fel44 


24d 


0.63/0.37 


0.44322(12) 


0.49968(13) 


0.40016(12) 


Pd/Fel45 


24d 


0.78/0.22 


0.49849(14) 


0.28713(13) 


0.24931(11) 


Pd/Fel46 


24d 


0.84/0.16 


0.34478(11) 


0.34920(11) 


0.03828(11) 


Pd/Fel47 


24d 


0.68/0.32 


0.49981(12) 


0.09886(11) 


0.05776(11) 


Pd/Fel48 


24d 


0.55/0.45 


0.15417(14) 


0.34369(13) 


0.15351(14) 


Pd/Fel49 


24d 


0.55/0.45 


0.40154(13) 


0.44069(12) 


0.00331(13) 



u 



site type 



A/aw 
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0.015 
0.132 
0.007 
0.004 
0.008 
0.160 
0.021 
0.014 
0.025 
0.025 
0.023 
0.034 
0.034 
0.067 
0.024 
0.018 
0.013 
0.019 
0.016 
0.012 
0.038 
0.023 
0.011 
0.001 
0.007 
0.012 
0.013 
0.009 
0.009 
0.006 
0.016 
0.007 
0.011 
0.009 
0.012 
0.005 
0.016 
0.010 
0.013 
0.009 
0.010 
0.009 
0.011 
0.010 
0.008 
0.013 
0.019 
0.008 
0.007 
0.009 
0.009 



B< 4) ,M< 5) 

B f\ M f> 
b£ 5) ,M< 5) 

b< 4) ,m< 6) 



B 



3 

(3) 

b< 3 \m< 6) 

Bf ,Mf 



B< 6) ,M< 6) 
B( 4 »,Mf» 

Bi 5) ,Mf 

M f>> 

B 3 7) ,M< 7) 



B 



(2) M (T) 



Bi 4) ,M« 

(6) iy/r(7) 



Bv ' ,M 



b| 4) ,m^ 

3(3) 



B^MV 

Br\Bi 7 

b{ 5 \b{ 6 



bKb 



b£ 5) ,B< 6 

b| 4) ,b| 7 

R (7) R (7 

b 3 ,a 3 
Bf,Bf 
B f\ B f 
Bi 4 \Bi 5 

b| 3) ,b| 6 

BfM 5 



Bf» 

Bf 

Bf» 

Bf» 

B< 4) 
B (i) 

°3 

B (4) 

°3 

B (6) 

°3 

B (D 

D 3 

B (6) 

D 3 

Bf» 

R(6) 



B 



3 

(2) 
3 

B (4) 
°3 

B f> 
B< 5) 



B3 



B 



B 3 4 
Bf 
Bf 

Bf 
Bf 
B (4 

°3 

Bf 
Bl 7 



B (4 

D 3 

Bf 
Bf 

Bf 
Bi 5 

B (6 
D 3 

B (7 

D 3 

(5 



(3 

m" 5 



(4 



M 



Mi 4 

Mf 
Mf 



M< 4 



M< 5 

M 5 

Mf 
Mf 
Mf 
Mf 
Mf 
Mf 
Mf 
Mf 
M* 1 



,M< B) > 
,M< 5) ) 



,Mf» 
,Mf» 
,M< 7) 



M ( 5 



((>) 



B3 
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,M< 3) 



0.003 
0.019 
0.001 
0.001 
0.003 
0.016 
0.001 
0.002 
0.002 
0.004 
0.003 
0.001 
0.003 
0.003 
0.000 
0.002 
0.002 
0.003 
0.002 
0.002 
0.003 
0.002 
0.003 
0.003 
0.003 
0.001 
0.004 
0.003 
0.003 
0.004 
0.002 
0.001 
0.005 
0.004 
0.005 
0.001 
0.003 
0.004 
0.003 
0.003 
0.003 
0.001 
0.002 
0.001 
0.002 
0.005 
0.002 
0.004 
0.004 
0.002 
0.004 
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atom 


Wyckoff 


S.O.f. 


x/ai a t 


y/a-iat 


z/ai a t 


Pd/Fel50 


24d 


0.57/0.43 0.44148(12) 


0.19141(12) 


0.09361(13) 


Pd/Fel51 


24d 


0.64/0.36 


0.40261(13) 


0.44054(13) 


0.18868(12) 


Pd/Fel52 


24d 


0.53/0.47 


0.43821(12) 


0.49800(13) 


0.09582(13) 


Pd/Fel53 


24d 


0.55/0.45 0.40714(14) 


0.44205(13) 


0.30777(13) 


Fe/Pdl54 


8c 


0.72/0.28 


0.15344(17) 


0.15344(17) 


0.15344(17) 


Fe/Pdl55 


24d 


0.58/0.42 


0.24809(13) 


0.30790(15) 


0.09506(16) 


Pd/Fel56 


24d 


0.64/0.36 0.30852(14) 


0.40455(14) 


0.25214(13) 


Pdl57 


24d 


0.66 


0.21000(14) 


0.24990(14) 


0.18636(15) 


Pdl57' 


24d 


0.34 


0.1911(3) 


0.2491(3) 


0.1549(3) 


Pd/Fcl58 


24d 


0.55/0.45 0.40290(15) 


0.25131(13) 


0.19141(15) 


Cr/A1159 


24d 


0.99/0.01 


0.0963(2) 


0.24942(19) 


0.0002(2) 


Cr/Fcl60 


24d 


0.60/0.40 


0.4035(3) 


0.2486(2) 


0.0013(2) 


Fe/Pdl61 


8c 


0.91/0.09 


0.3445(3) 


0.3445(3) 


0.3445(3) 


Fe/Crl62 


24d 


0.80/0.20 


0.3470(2) 


0.3444(2) 


0.1534(2) 


Fe/Crl63 


24d 


0.70/0.30 0.24888(18) 


0.1904(2) 


0.0955(2) 


All 64 


24d 


1.00 


0.0602(4) 


0.2148(4) 


0.0364(4) 


A1165 


24d 


1.00 


0.4990(4) 


0.3805(4) 


0.3086(4) 


A1166 


24d 


1.00 


0.2270(3) 


0.4421(4) 


0.0947(4) 


A1167 


24d 


1.00 


0.1915(4) 


0.5005(4) 


0.1184(4) 


A1168 


24d 


1.00 


0.2488(4) 


0.4763(4) 


0.1540(4) 


A1169 


24d 


1.00 


0.4037(4) 


0.3463(4) 


0.1326(4) 


All 70 


8c 


1.00 


0.1900(4) 


0.1900(4) 


0.1900(4) 


A1171 


24d 


1.00 


0.2497(4) 


0.2140(5) 


0.1531(5) 


All 72 


8c 


1.00 


0.3092(5) 


0.3092(5) 


0.3092(5) 


All 73 


24d 


0.53 


0.2086(8) 


0.1610(8) 


0.1295(8) 


A1173' 


24d 


0.47 


0.190(1) 


0.1909(9) 


0.1204(9) 


All 74 


24d 


0.72 


0.0964(5) 


0.2698(5) 


0.0594(5) 


All 74' 


24d 


0.54 


0.0515(6) 


0.2840(6) 


0.0041(6) 


All 75 


24d 


0.61 


0.2834(6) 


0.5380(6) 


0.1334(6) 


A1175' 


24d 


0.39 


0.3139(9) 


0.4941(9) 


0.1078(9) 


All 76 


24d 


0.69 


0.2869(7) 


0.1605(7) 


0.1314(7) 


A1176' 


24d 


0.31 


0.1151(15) 


0.3102(16) 


0.1911(16) 


All 77 


24d 


0.45 


0.2306(11) 


0.250(1) 


0.0973(12) 


All 77' 


24d 


0.55 


0.2659(9) 


0.2487(8) 


0.0952(9) 


All 78 


24d 


0.52 


0.3678(11) 


0.2852(11) 


0.1503(11) 


A1178' 


24d 


0.48 


0.3790(11) 


0.3079(12) 


0.1827(11) 


All 79 


24d 


0.64 


0.0964(7) 


0.3458(7) 


0.3218(6) 


A1179' 


24d 


0.36 


0.0996(12) 


0.3508(12) 


0.3774(11) 


A1180 


24d 


0.62 


0.5259(8) 


0.3437(8) 


0.2504(7) 


A1180' 


24d 


0.38 


0.3484(12) 


0.2503(11) 


-0.0147(12) 


A1181 


24d 


0.48 


0.3824(12) 


0.3813(12) 


0.1876(12) 


A1181' 


24d 


0.52 


0.3683(11) 


0.4001(11) 


0.1561(11) 


A1182 


24d 


0.57 


0.347(1) 


0.367(1) 


0.286(1) 


A1182' 


24d 


0.43 


0.3129(13) 


0.3786(13) 


0.3085(13) 


A1183 


24d 


0.23 


0.5371(17) 


0.2480(16) 


0.0218(17) 


A1183' 


24d 


0.51 


0.5538(7) 


0.2776(7) 


0.0321(7) 


A1184 


24d 


0.76 


0.4253(4) 


0.0003(5) 


0.0002(5) 


A1184' 


24d 


0.24 


-0.0014(14) 


0.4007(12) 


-0.0007(15) 


A1185 


24d 


0.66 


0.4365(9) 


0.2128(9) 


0.0366(9) 


A1185' 


24d 


0.34 


0.3913(17) 


0.2214(17) 


0.0618(17) 


A1186 


24d 


0.61 


0.4403(6) 


0.4773(5) 


0.0363(6) 



u 



site type 



A/aw 
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0.005 
0.014 
0.010 
0.012 
0.011 
0.009 
0.015 
0.006 
0.009 
0.015 
0.008 
0.019 
0.057 
0.021 
0.002 
0.009 
0.016 
0.010 
0.013 
0.008 
0.014 
0.000 
0.022 
0.029 
0.029 
0.029 
0.011 
0.011 
0.017 
0.017 
0.038 
0.038 
0.036 
0.036 
0.057 
0.057 
0.029 
0.029 
0.031 
0.031 
0.059 
0.059 
0.058 
0.058 
0.015 
0.015 
0.009 
0.009 
0.067 
0.067 
0.007 



,M< 7) 



B Wb< 8 \m< 4 \m< 6 >> 

p4 4 \p4 7 \m( 3) 

(4) 

3 

(1) 



Bl 2) ,B| 4) ,Ni 



B^M^M^MD 
B< 3) ,B< 8) M 5) ,M( 6) > 

b| 3 \b| 4) ,Iv4 7 \m< 7) > 

B< 3) ,M< 5) > ' 

b| 3) ,m< 5) } 
b< 3 \b< 6) ,m< 6) 



,Mf») 
,M« 7) > 

f(5)\ 



fC)\ 



m< 3 \m< 4 \m< 6 > 
b| 7) ,m< 4) 
b< 2) ,m< 7) 



Ml 

) 

m( 6) ,m< 7) ) 

00 \ 



b< 3) ,m< 3 \m< 6) 



b' 4) 



M 



3 

Mf 
M< 3 

Mf 

Mf 
M< 5 
M< 5 
Mf 
Mf 
Mf 



Mr' 



(s) ,Mf ,M< 6) > 
,M^ 5) > 
,M< 7) ) 
> 
> 

M 7) > 

M^) 
,M< 7) ) 



Mf 
Mf 
M (3 



M 



2 

Mf 

Mf 

Mf 
Mf 

(7 

Mf 

(3 

Mf 
Mf 



Ml, 



Ml, 



Mf 
Mf 
Mf 
M (r 

(7 
2 

Mf 
Mf 
Mf 
Mf 

(6 

Mf 
Mf 



Ml, 



Ml 



Mf 
Mf 
Mf 
Mf 
Mf 
Mf 
Mf 
Mf- 



0.002 
0.003 
0.003 
0.003 
0.002 
0.002 
0.002 
0.006 
0.042 
0.002 
0.001 
0.002 
0.002 
0.002 
0.001 
0.002 
0.002 
0.001 
0.001 
0.002 
0.001 
0.002 
0.002 
0.000 
0.009 
0.003 
0.003 
0.033 
0.004 
0.013 
0.006 
0.003 
0.004 
0.007 
0.004 
0.009 
0.001 
0.012 
0.004 
0.008 
0.003 
0.005 
0.002 
0.005 
0.002 
0.011 
0.002 
0.019 
0.005 
0.015 
0.001 
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atom wyckoff s.o.f. x/a lat y/aiat z/a lat U site type A/g; at 

(continued from the previous page) 
A1186' 24d 0.39 0.4193(9) 0.4689(8) 0.0510(9) 0.007 (M^) 0.018 

A1187 24d 0.59 0.4658(9) 0.4381(9) 0.2866(9) 0.050 (M< 7) ) 0.004 

A1187' 24d 0.41 0.4423(13) 0.4010(13) 0.2735(12) 0.050 (M< 4) ) 0.004 

(end of table) 
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